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ABSTRACT 


Absorbance  and  steady-state  and  time-resolved  fluorescence  measurements  were 
employed  to  examine  the  mechanism(s)  of  excited  singlet  state  quenching  by  nitroxides  in  a 
series  of  nitroxide-fluorophore  adducts.  This  work  establishes  the  following:  1)  the  absorption 
and  emission  energies  of  the  fluorophores  are  unaffected  by  the  presence  of  the  nitroxide 
substituent(s),  and  the  residual  emission  that  is  observed  from  the  adducts  arises  from  the 
locally  excited  singlet  of  the  fluorophore,  not  from  charge  recombination;  2)  rate  constants  for 
intramolecular  quenching  by  the  nitroxides  (kq)  are  high  (108  -1010  s*1)  and  decrease 
significantly  with  increasing  nitroxide  to  fluorophore  distance  -  however,  relatively  high  rates  of 
quenching  (MO8  s"1)  are  observed  over  distances  as  great  as  12  A;  3)  Forster  energy  transfer 
does  not  contribute  significantly  to  the  quenching  due  to  the  low  values  for  the  spectral  overlap 
integrals;  4)  the  kq's  do  not  increase  proportionally  to  the  solvent-dependent  increases  in  the 
Dexter  overlap  integral,  indicating  that  energy  transfer  by  the  Dexter  mechanism  is  not 
responsible  for  the  quenching;  5)  the  values  of  kq  show  no  obvious  correlation  with  the 
calculated  free  energies  for  photo-induced  electron  transfer,  suggesting  that  this  quenching 
pathway  is  also  unimportant;  6)  for  hematopotphyrin-nitroxide  adducts,  which  contain  a 
fluorophore  whose  singlet  energy  is  below  that  of  the  first  excited  state  energy  of  the  nitroxide 
(thus  precluding  energy  transfer),  significant  rates  of  quenching  are  still  observed;  7)  for 
compounds  with  similar  nitroxide-fluorophore  distances,  an  approximately  linear  correlation  is 
observed  between  the  kq’s  of  the  paramagnetic  compounds  and  the  non-radiative  rate  constants 
of  the  diamagnetic  reference  compounds,  suggesting  that  the  nitroxide  moiety  catalyses  a  pre¬ 
existing  non-radiative  pathway  in  the  fluorophore.  These  results  indicate  that  the  quenching 
arises  through  electron  exchange  which  causes  relaxation  of  the  (local)  singlet  state  to  the  triplet 
and/or  ground  state  of  the  fluorophore. 


Introduction 


Diffusional  quenching  by  stable  nitroxyl  radicals  of  excited  singlet1,  doublet 2  and 
triplet1^*3'10  states  and  of  excimers  11  has  been  studied  extensively  over  the  past  twenty  years. 
The  goal  of  much  of  this  work  has  been  to  understand  the  mechanism(s)  through  which  excited 
states  are  quenched  by  paramagnetic  species.  Recently  however,  fluorescence  quenching  by 
nitroxides  has  also  become  an  important  tool  with  which  to  probe  the  structural  and  dynamical 
properties  of  membranes12  and  micelles.13 

While  it  has  been  generally  accepted  that  nitroxide  quenching  of  singlet  states  can  result 
from  an  electron  exchange-induced14  intersystem  crossing  to  the  triplet lb>e  or  internal 
conversion  to  the  ground  state,  le'8  significant  contributions  to  the  quenching  may  also  arise 
from  chargelf  and/or  energy  transfer. la»12d  Green  et  al.lb  argued  that  charge  (electron)  transfer 
is  unimportant  owing  to  a  lack  of  solvent  dependence  of  the  quenching  rate  constant.  However, 
as  noted  by  Chattopadhyay  et  al.,lf  the  quenching  occurs  at  the  diffusion  limit  and  thus  may  not 
be  differentiated  by  solvent  polarity.15  Green  et  al.lb  also  suggested  that  Forster  energy 
transfer16  is  insignificant  because  of  the  low  extinction  coefficients  of  the  nitroxides  in  the  400- 
550  nm  spectral  region.  Recent  work  by  Puskin  et  al.12d  suggests  that  energy  transfer  can  be 
important  and,  in  some  cases,  can  extend  the  quenching  radius  to  as  much  as  10  A, 
approximately  twice  the  interaction  distance  originally  calculated  by  Green  et  al.lb  Furthermore, 
energy  transfer  by  electron  exchange17  cannot  be  excluded  a  priori  for  compounds  having 
singlet  energies  greater  than  the  lowest  excited  state  energy  of  the  nitroxide  (~2.6  eV  for  the 
piperidinyl  nitroxides). 

Until  recently,  information  on  the  efficiency  of  intramolecular  quenching  of  excited 
states  by  nitroxides  was  unavailable.  We  have  shown  that  intramolecular  quenching  of  excited 
singlets  is  highly  efficient,18  and  that  covalently  linked,  nitroxide-fluorophore  adducts  can  be 
employed  as  very  sensitive  optical  probes  of  radical/redox  reactions.18'21  These  adducts  also 
offer  a  means  of  studying  the  mechanism(s)  of  nitroxide  quenching  in  more  detail  and  of 
establishing  the  distance  dependence  and  possible  importance  of  "through-bond"15  vs  "through- 


space"  interactions  in  the  quenching.  This  information  is  essential  for  the  construction  of  better 
optical  radical  sensors,18*20  for  a  more  complete  analysis  of  nitroxide  quenching  in  biophysical 
studies12*13  and  for  an  understanding  of  the  possible  influence  of  stable  or  transient  radicals  on 
excited  states  in  biological  and  chemical  systems. 

In  this  study,  we  examine  the  photophysical  properties  of  a  number  of  nitroxide- 
fluorophore  adducts  and  their  diamagnetic  analogues.  By  employing  compounds  having 
fluorophores  with  differing  singlet  energies  and  a  range  of  nitroxLe-fluorophore  distances,  we 
attempt  to  assess  the  relative  contributions  of  energy  transfer,  electron  transfer,  and  electron 
exchange-induced  relaxation  to  the  intramolecular  quenching  of  singlet  states  by  nitroxides. 

Results  and  Discussion 

A.  Absorption.  Excitation,  and  Emission  Spectra 

Table  1  summarizes  the  absorption  and  emission  energies  of  I-IV  (Fig.  1)  in  various 
solvents.  In  water,  methanol,  acetonitrile,  dioxane  and  hexane,  the  absorption  spectra  of  the 
nitroxide-fluorophore  adducts  were  indistinguishable  from  the  sum  of  the  individual  spectral 
contributions  of  the  fluorophore  and  nitroxide  substituents.  Other  than  the  loss  of  the  weak 
nitroxide  contribution  in  the  visible  (430-480  nm,  £=10-13  M'lcm'l)  and  ultraviolet  (-246  nm, 
e=2xl03  M*  1cm*1),  the  spectra  of  the  O-acetyl  and  methyl  ester  derivatives  (Fig.  1)  were 
identical  to  those  of  the  corresponding  nitroxides.  No  evidence  for  optical  charge  transfer  bands 
or  enhanced  ground  to  triplet  absorption14  by  the  fluorophore  could  be  found  in  the  spectra  of 
the  nitroxide  adducts.  While  the  absorption  spectra  of  the  fluorophores  were  largely  insensitive 
to  solvent  polarity,  the  lowest  energy  band  (n-7t*)  of  the  nitroxides22  shifted  substantially  to  the 
blue  with  increasing  solvent  polarity  as  previously  noted.23  The  influence  of  this  shift  on 
singlet  state  quenching  will  be  addressed  below  (Section  D). 

The  fluorescence  excitation  spectra  of  the  nitroxide  adducts  matched  the  absorption 
spectra  of  the  fluorophore  substituents  and  were  indistinguishable  from  the  excitation  spectra  of 
the  corresponding  diamagnetic  derivatives. 
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The  lineshapes  of  the  emission  spectra  were  identical  for  the  nitroxide  adducts  and  the 
corresponding  diamagnetic  derivatives,  and  generally  showed  a  mirror-image  relationship  to  the 
excitation  spectra  (Fig.  2).  With  increasing  solvent  polarity,  identical  red-shifts  in  the  emission 
maxima  were  observed  for  the  nitroxide  and  diamagnetic  adducts  of  a  given  substitution.  These 
results  indicate  that  the  residual  emission  from  the  nitroxide  adducts  is  arising  from  the  same 
locally  excited  singlet  state  of  the  fluorophore  substituent  as  in  the  diamagnetic  derivatives  and  is 
not  produced  through  charge  recombination.24  No  evidence  for  exciplex  emission  from  the 
paramagnetic  adducts  could  be  found  in  either  polar  (acetonitrile)  or  nonpolar  (hexane)  solvents. 

B.  Quantum  Yield  and  Lifetime  Measurements 

While  the  presence  of  the  nitroxyl  radical  substituent(s)  did  not  influence  the 
absorption  and  emission  energies  of  the  fluorophores  in  the  adducts,  the  fluorescence  quantum 
yields  and  lifetimes  of  the  nitroxide  adducts  were  reduced  substantially  relative  to  their 
corresponding  diamagnetic  derivatives  in  both  polar  and  nonpolar  solvents  (Table  2,  Fig.  3). 
The  fluorescence  decay  of  the  nitroxide  adducts  could  be  well  fit  to  a  sum  of  two  exponentials 
but  was  dominated  by  the  short-lived  component  (Fig.  3).  The  small,  long-lived  component 
usually  represented  <5%  of  the  total  amplitude.  We  attribute  this  component  to  minor 
diamagnetic  contaminants,  based  on  assessments  of  compound  purity. 

With  few  exceptions,  the  fluorescence  decay  of  the  diamagnetic  derivatives  could  be  fit 
as  a  single  exponential  process  (Fig.  3).  The  O-acetyl  and  methyl  ester  derivatives  of  a  given 
substitution  showed  very  similar  lifetimes  and  quantum  yields  (Table  2).  These  results  indicate 
that  the  rate  constants  for  radiative  and  non-radiative  decay  of  the  singlet  do  not  differ  signifi¬ 
cantly  between  the  O-acetyl  and  the  methyl  ester  derivatives  (Fig.  1).  This  leads  us  to  conclude 
that,  unlike  the  situation  for  aliphatic  and  aromatic  amines,25  intramolecular  quenching  (via 
electron  transfer)  by  the  hindered,  O-substituted  hydroxylamines  is  unimportant.  Preliminary 
work  on  other  O-substituted  hydroxylamines  is  consistent  with  this  conclusion.19  In  contrast, 
the  unsubstituted  hydroxylamine  does  appear  to  quench  partially.19-26 


The  steady-state  and  time-resolved  fluorescence  measurements  of  the  paramagnetic  and 


diamagnetic  forms  can  be  related  by 


(1) 

(hq  =  **  , 

kf  +  knr  =  ^ 

(2) 

^  kr  +  knr  +  kq 

kr  +  krr  +  kq=:-L 

where  <j)d  (<J>p)  and  t<]  (tp)  are  the  quantum  yields  and  lifetimes  of  the  diamagnetic  (paramagnetic) 
compounds,  respectively,  kr  and  knr  are  the  rate  constants  for  radiative  and  non-radiative  decay 
of  the  diamagnetic  compounds,  and  kq  is  the  intramolecular  (non-radiative)  quenching  rate 
constant  attributable  to  the  nitroxide  moiety.  If  kr  and  knr  are  unaltered  by  the  presence  of  the 
nitroxide,  then 


(3) 

(4) 


<t>d  _T<1 
0p  TP 


k^-L-JL 

^  'Cp  Td 


The  measured  ratios  Xd/xp  and  <J>dAJ>p  are  equal  within  the  uncertainties  (Table  2), 
indicating  that  this  condition  applies.  Errors  are  largest  in  the  determination  of  <|>p  since  traces  of 
highly  fluorescent  impurities  can  produce  significant  overestimates  of  such  low  quantum  yields, 
leading  to  low  <bd/^>p  ratios.  For  this  reason,  rate  constants  for  intramolecular  quenching,  kq 
(Table  3)  were  calculated  from  the  lifetime  data  employing  Eqn.  4,  except  in  the  case  of 
compound  IA  where  lifetimes  were  too  short  to  be  resolved  by  the  instrument.  For  this 
compound  the  expected  lifetimes  were  computed  from  Eqn.  3  using  the  quantum  yield  data  and 
the  lifetime  of  the  diamagnetic  O-acetyl  (IB).  Lifetime  data  show  <3%  of  a  long-lived 
component  for  IA.27 

In  aprotic  solvents  such  as  acetonitrile  and  hexane,  kq  increased  from  “108  s*1  to 
10 10  s*1  in  the  order  ID  <  IIA  <  IDA  <  IVA  <  IA.  This  sequence  reflects  not  only  the 


differences  in  nitroxide-fluorophore  distance  but  also  the  type,  number  and  position  of  the 
naphthalene-nitroxide  linkage(s)  (vide  infra).  With  the  exception  of  compounds  IA  and  ID,  no 
significant  differences  were  observed  between  kq's  measured  in  polar  (acetonitrile)  and 
nonpolar  (hexane)  solvents.  Surprisingly,  the  kq’s  for  compounds  IA  and  ED  actually  increased 
slightly  with  decreasing  solvent  polarity  (Table  3). 

C.  Dependence  of  kq_on  Distance.  Linkage  Position  and  Type 

Because  the  distance  and  orientation  of  the  nitroxide  with  respect  to  the  fluorophore  are 
not  rigidly  fixed  in  these  compounds,  the  molecular  mechanics  program  AMBER28  was 
employed  to  search  for  the  conformations  of  IA,  IIA  and  MA  most  likely  to  exist  in  solution. 
The  strategy  for  this  search  is  described  in  Experimental. 

IA  and  IIIA  were  each  found  to  have  three  unique  minimal  energy  conformations 
(Table  4).  However,  the  average  energy  of  the  low  energy  conformations  of  IIIA  was  5.6 
kcal/mol  lower  than  that  of  IA.  The  higher  overall  energy  of  IA  can  be  attributed  to  a  sterically- 
induced  distortion  of  the  naphthalene  ring  and  subsequent  effects  on  the  dihedral  angles  within 
the  ring  and  between  the  ring  and  the  carboxylate  group.  This  distortion  was  a  result  of 
unfavorable  steric  interactions  between  the  carbonyl  oxygen  of  the  carboxyl  group  and  the 
hydrogen  at  position  8  on  the  naphthalene  ring.  A  rotation  of  the  naphthalene  ring  by  180* 
relative  to  the  carboxylate  group  resulted  in  greater  steric  interactions  and  a  more  distorted 
structure.  Neither  of  the  other  compounds  studied  displayed  this  distortion. 

In  contrast  to  IA  and  IDA,  IIA  exhibited  8  low  energy  conformations  (Table  4).  A 
single  conformer  exhibited  both  the  lowest  energy  and  shortest  ring  to  nitrogen  distance.  Seven 
other  conformers  had  an  average  energy  that  was  only  0.7  kcal/mol  higher.  Because  these 
conformers  constituted  a  large  set  of  thermally  accessible  states,  we  were  forced  to  include  them 
in  the  calculation  of  average  distances.  The  energies  and  center-to-center  distances  of  the 
energetically  significant  conformations  of  IA,  IIA,  and  IIIA  are  summarized  in  Table  4. 
Coordinates  for  each  of  these  conformations  are  provided  in  the  supplementary  material. 


The  fact  that  the  short-lived  fluorescence  component  of  each  of  these  three  compounds 
could  be  well  fit  as  a  single  exponential  decay  suggests  either  that  these  conformers  interconvert 
rapidly  on  the  time  scale  of  the  quenching  (ns  to  sub-  ns),  or  that  conformations  populated 
during  the  lifetime  of  the  singlet  do  not  differ  significantly  in  their  ability  to  quench  the  singlet 
state. 

The  Forster  orientation  factor,  k2,  was  calculated  for  the  conformers  of  IA-IIIA  (Table 
4)  as  a  fairly  simple  means  to  compare  the  relative  orientations  of  naphthalene  and  nitroxide 
groups.  The  absence  of  large  variations  in  k2  suggests  that  differences  in  the  relative 
orientation  of  the  two  groups  is  probably  not  the  primary  factor  controlling  the  variations  in 
cuenching  rates  among  these  compounds. 

The  rates  of  energy  transfer  (Forster  or  Dexter)  and  electron  transfer  have  well-defined 
distance  dependencies.  While  the  rate  of  Forster  energy  transfer  follows  an  R*6  distance 
dependence,16-29  the  rates  of  Dexter  transfer  and  electron  transfer  are  expected  to  decline 
exponentially  with  distance.15-17*30*31  However,  for  aprotic  solvents,  the  observed  quenching 
rates  bear  no  simple  relationship  to  the  nitroxide-naphthalene  distances  (Tables  3  and  4),  as  is 
illustrated  by  the  following  points:  1)  while  a  comparison  between  IA  and  ID  suggests  that 
increasing  R  by  4  A  causes  kq  to  decrease  by  about  two  orders  of  magnitude,  simply  shifting 
the  site  of  substitution  from  the  a  (IA)  to  the  |3  (HLA)  position  on  the  naphthalene  ring  produces 
an  -20-fold  decrease  in  kq,  although  the  values  of  R  for  IA  and  mA  differ  by  only  0.76  A;  2) 
insertion  of  a  methylene  group  into  the  linkage  at  the  a  position  (HA)  decreases  kq  by  -40-fold, 
although  the  greater  flexibility  of  this  linkage  allows  a  shorter  R  than  either  IA  or  HLA.  As  we 
show  below,  the  rather  large  variations  in  kq  observed  for  compounds  having  similar  nitroxide- 
naphthalene  distances  (e.g.  IA-IVA)  cannot  be  explained  by  the  differences  in  the  singlet  energy 
of  the  naphthalenes  or  by  variations  in  the  free  energy  for  photo-induced  electron  transfer. 
Instead  these  variations  appear  to  result  from  a  relationship  between  kq  and  the  non-radiative  rate 
constant  of  the  parent  (diamagnetic)  compound.  The  presentation  and  interpretation  of  this 
relationship  is  provided  in  section  D  3. 


D.  Mechanisms  of  Quenching 


Several  mechanisms  could  contribute  to  singlet  state  quenching  in  these  nitroxide 
adducts.  We  address  here  the  possible  contributions  of  Forster  (dipole- dipole)  and  Dexter 
(electron  exchange)  energy  transfer,  electron  transfer,  and  electron  exchange-induced  relaxation 
to  the  quenching.  In  principle,  kq  reflects  the  sum  of  rate  constants  for  each  of  these  potential 
relaxation  pathways: 

(5)  kq  =  1<FT  +  koT  +  key  +  kex 


1.  Energy  Transfer 

In  the  light  of  conflicting  reportsla’lb>lf*12d  as  to  the  importance  of  Forster  and 
Dexter  transfer  in  the  quenching  of  singlet  states  by  nitroxides,  we  have  re-examined  this 
question.  Since  no  emission  is  observed  from  the  first  excited  state  of  the  nitroxide,  direct 
spectroscopic  evidence  for  energy  transfer  was  not  obtainable.  Instead,  we  compared  the 
observed  quenching  rates  of  IA,  HA,  and  HIA  to  the  Forster  (Jp)  and  Dexter  (Jd)  overlap 
integrals  calculated  from  the  fluorescence  spectra  of  EB,  KB  and  KIB,  and  the  absorbance 
spectrum  of  4- hydroxy-2,2, 6,6-tetramethyl-l-piperdinyloxy  (TEMPOL): 
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where  F(X)  is  the  corrected  fluorescence  intensity  of  the  fluorophore  and  £a(X)  is  the  extinction 
coefficient  of  the  acceptor  (TEMPOL)  at  wavelength  X. 

The  weak  nitroxide  absorption  band  overlaps  with  the  red  edge  of  the  singlet  emission 
of  the  naphthyl  group.  This  absorption  is  solvent  dependent23  (Fig.  4),  shifting  to  the  blue 
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with  increasing  solvent  polarity  and  thus  increasing  the  values  of  Jf  and  Jd  •  The  degree  of 
spectral  overlap  is  further  enhanced  in  polar  solvents  for  compounds  IA  and  HLA  whose 
emissions  shift  to  the  red  with  increasing  solvent  polarity  (Table  1). 

We  have  used  these  solvent-dependent  changes  in  the  overlap  integrals  to  assess 
whether  the  quenching  rates  vary  in  the  manner  expected  for  energy  transfer.  The  rate  constant 
for  Forster  transfer  is  given  by, 29 

(7a)  kFT  =  JfK2|i>d  8,71  x  io23 

R^rvkd 


where  k2  is  the  orientation  factor  (Table  4),  <j>d  and  xp  are  the  quantum  yield  and  lifetime, 
respectively,  of  the  appropriate  diamagnetic  O- acetyl  derivative  (Table  3),  R  is  the  donor- 
acceptor  distance  (Table  4)  ind  T|  is  the  solvent  refractive  index. 

The  rate  constant  for  Dexter  transfer  is  given  by,15b 


where  H  is  the  coupling  matrix  element,  Ho  is  the  pre-exponential  factor  and  L  is  the  Bohr 
radius. 

Both  of  these  theories  predict  that  for  constant  donor-acceptor  distance(s)  (and 
orientation(s))  the  rate  constants  for  transfer  will  increase  linearly  with  increasing  values  of  the 
spectral  overlap  integrals.  However,  log  plots  of  kq  vs  Jf  and  Jd  for  IA,  IIA,  and  IIIA  (Fig.  5) 
clearly  do  not  show  this  relationship  (with  the  possible  exception  of  IDA  in  the  Forster  plot 
(Fig.  5)).  Furthermore,  equation  7a  can  be  employed  to  calculate  kFT  from  R  and  k2  values 
(Table  4);  these  predicted  Forster  transfer  rates  (Table  5)  are  more  than  an  order  of  magnitude 
lower  than  the  measured  values  of  kq.  Moreover,  choosing  a  value  of  2/3  for  k2, 
corresponding  to  randomized  donor-acceptor  orientations,29  increases  Icft  to  a  maximum  of 
only  -15%  of  the  observed  kq  values.  Even  in  the  case  of  anthracene,  where  the  spectral 


overlap  with  TEMPOL  is  maximal  for  the  family  of  linear  aromatic  ring  systems,  we  calculate 
that  the  quenching  radius  for  Forster  energy  transfer  extends  to  only  9  A  (in  ethanol). 

We  conclude  that  energy  transfer  through  either  a  Forster  or  Dexter  mechanism  cannot 
account  for  the  observed  magnitude  and  (lack  of)  variation  in  the  intramolecular  quenching  rate 
constants. 

2.  Electron  Transfer 

For  electron  transfer  to  be  operative,  the  energy  of  the  excited  singlet  must  be  sufficient 
to  drive  the  relevant  redox  couple.30*31  In  the  absence  of  Coulombic  interactions,  the 
thermodynamic  driving  force,  AG«„  is  given  by30 

(8)  AGoo  —  -Eqo  +  Eox  "  Ered 

where  Eoo  is  the  00-excitation  energy  of  the  fluorophore,  and  Eox  and  Ered  are  the  oxidation 
potential  of  the  electron  donor  and  the  reduction  potential  of  the  acceptor,  respectively. 

Assuming  a  simple  continuum  model  for  solute-solvent  interaction,  additional  terms 
for  the  Coulombic  attraction  between  the  separated  charges,  AGCoul.  and  for  the  free  energy  of 
solvation  of  the  ion  pair,  AGsolv  must  be  included.  The  total  change  in  free  energy  is  then 
given  by:30*31 


(9)  AGet  =  AG*  +  AGcoul  +  AGsolv 


(10a)  AGcoa  =  -e2/eR 

(10b)  AGsov  =  ^(J-  +  J-)(— 1— -1) 

2  vd  fa/leMeCN  el 


and  e  is  the  electron  charge,  £  is  the  solvent  dielectric  constant,  £MeCN  is  the  dielectric  constant 
of  acetonitrile  in  which  the  redox  potentials  were  measured,  R  is  the  distance  between  donor  and 
acceptor  (Table  4),  r<i  and  tq  are  donor  and  acceptor  radii  espectively.  Coulombic  and  solvation 
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adjustments  to  AG„  are  greatest  (up  to  +0.8  eV)  in  hexane  (e  =  1.88)  and  increase  steeply  with 
distance  in  this  solvent.  For  the  polar  solvents  acetonitrile  and  water,  AGcoul +  AGsolv  < 

0.1  eV  and  exhibits  little  distance  dependence. 

For  IA-IVA,  the  naphthalene  moiety  may  act  as  either  an  electron  acceptor  or  donor 
resulting  in  oxidation  or  reduction,  respectively,  of  the  nitroxide;  these  reactions  are  shown 
diagramatically  in  eqns.  11  and  12: 

no’  (ID 

NO  (12) 

acetonitrile  both  reduction  and  oxidation  of  the  nitroxide  are  thermodynamically  feasible  but  that 
the  oxidative  pathway  (rxn  12)  is  preferred  in  all  cases  by  >0.15  eV.  In  hexane,  the  lack  of 
solvent  stabilization  makes  reduction  of  the  nitroxide  (rxn  11)  energetically  unfavorable  and 
reduces  the  driving  force  for  the  oxidative  pathway  (AGet  >  -0.4  eV).  In  aqueous  solution, 
nitroxide  reduction  (rxn  1 1)  has  the  greater  driving  force.  Although  these  calculations  suggest 
that  the  direction  of  electron  transfer  would  reverse  in  water,  the  quenching  rate  constants  do  not 
differ  greatly.  If  electron  transfer  is  the  dominant  pathway  for  quenching,  transfer  in  both 
directions  must  be  occurring  with  equal  facility. 

Within  individual  series  of  compounds  our  measured  values  of  kq  do  not  show  any 
obvious  correlation  with  AGet-  For  example,  the  values  of  kq  for  ILA  are  equal  in  acetonitrile 
and  hexane,  although  electron  transfer  is  calculated  to  be  exergonic  in  the  first  solvent  and 
endergonic  in  the  latter.  Other  discrepancies  are  also  apparent  In  particular,  the  kq's  for  IUA 
and  IV A  do  not  vary  significantly  between  acetonitrile  and  hexane  despite  the  large  (-0.8  eV) 
differences  in  AGet  in  these  two  solvents.  Additionally,  in  polar  solvents  quenching  of  niA, 
(^-substituted)  is  slower  than  LA  (a-substituted)  by  an  order  of  magnitude  although  its  driving 
force  is  greater. 


Overall,  these  data  are  inconsistent  with  the  notion  that  electron  transfer  acts  as  the 
dominant  pathway  for  singlet  quenching I5aj0,3l#  The  lack  of  solvent  dependence  does  not,  by 
itself,  rule  out  electron  transfer  since  the  rate  of  transfer  is  a  function  of  the  activation  energy, 
AG$,  which  includes  an  additional  solvent  dependence  in  the  form  of  the  solvent  reorganization 
energy,  XS.30>31  Oevering  et  al.^03  have  demonstrated  that  the  effect  of  solvent  reorganization 
on  AG$  may  in  some  cases  largely  compensate  for  the  solvent  dependence  included  in  AGej. 
We  were  unable  to  calculate  AG$  precisely  for  these  compounds  owing  to  uncertainties  in 
sensitive  parameters  such  as  the  internal  (nuclear)  reorganization  energy  and  the  effective  radius 
of  the  nitroxide  group.  Nevertheless,  no  reasonable  values  for  these  parameters  lead  to  a 
prediction  of  faster  rates  in  hexane  than  in  acetonitrile  as  we  have  observed  for  several 
compounds  (IA,D).  Based  on  these  observations,  we  tentatively  conclude  that  electron  transfer 
does  not  play  a  major  role  in  the  quenching.  The  complete  lack  of  radical  ion  formation  during 
singlet  state  quenching  by  nitroxides  in  polar  solvents  lf*8  provides  additional  support  for  this 
view.  However,  owing  to  the  uncertainties  in  the  simple  continuum  model  and  the  redox 
potentials  for  the  nitroxides,2»3234  3$  yet>  we  are  hesitant  to  exclude  unequivocally  a 
contribution  of  electron  transfer  to  the  quenching. 


3.  Electron  Exchange 

Electron  exchange  between  the  naphthoate  singlet  and  nitroxide  doublet  can  produce 
rapid  relaxation  of  the  excited  state.1*14  Available  relaxation  routes  are  as  follows:1  (Dexter) 
energy  transfer  (13a),  intersystem  crossing  to  triplet  (13b),  or  internal  conversion  (13c) 
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As  discussed  above  Dexter  energy  transfer  is  an  unlikely  relaxation  pathway. 

However,  to  eliminate  energy  transfer  from  consideration  and  focus  on  exchange-induced 
intersystem  crossing  (13b)  or  internal  conversion  (13c),  preliminary  data  were  collected  on 
several  hematoporphyin-nitroxide  adducts  (Fig.  6).  For  these  compounds,  the  energy  of  the 
porphyin  singlet  (~2.0  eV)  is  well  below  that  of  the  pyrrolidinyl  nitroxide  doublet  (~2.9  eV) 
thus  eliminating  the  possibility  of  Dexter  (or  Forster)  energy  transfer.  AGet  for  these  adducts  is 
also  calculated  to  be  small  (-0.18  eV,  nitroxide  oxidation)36.  Nevertheless,  in  methanol,  the 
rate  constants  for  quenching  in  the  A  isomer,  6.6x10s  s*1,  and  the  B  isomer,  7.2x10s  s*1,  are 
comparable  to  those  of  HA.  These  results  suggest  that  an  electron-exchange  induced 
intersystem  crossing  to  the  triplet  or  internal  conversion  to  the  ground  state  is  occurring. 

Further  support  for  an  exchange  induced  "local"  relaxation  is  provided  by  the 
correlation  presented  in  Fig.  7.  Here  we  have  plotted  the  kq's  of  the  paramagnetic  compounds 
against  the  non-radiative  rate  constants  of  their  corresponding  diamagnetic  O-acetyl  and  methyl- 
ester  derivatives.  With  the  exception  of  IVA  in  methanol,  compounds  having  similar  nitroxide- 
fluorophore  distances  (e.g.  IA,  IIA-IVD)  exhibit  an  approximately  linear  relationship  between 
these  parameters.  This  correlation  is  interpreted  to  result  from  a  fixed  enhancement  (~55  -  fold) 
by  the  nitroxide  of  a  pre-existing,  non-radiative  relaxation  pathway  in  the  fluorophore.  Because 
non-radiative  decay  in  most  aromatic  hydrocarbons  proceeds  almost  exclusively  through 
intersystem  crossing  (Ermoiev's  rule14c),  this  enhancement  probably  reflects  a  lifting  of  the 
spin-restriction  to  singlet-triplet  transition  brought  about  through  the  coupling  of  the  doublet  and 
excited  singlet  states.14  We  suspect  that  the  relatively  minor  deviations  from  this  correlation 
result  from  variations  in  the  degree  of  coupling  among  these  derivatives  due  to  small  differences 
in  distance,  orientation,  and  structure  (Fig.  1,  6). 

If  intersystem  crossing  (13b)  competes  effectively  with  triplet  quenching  by  the 
nitroxide  (13b  — >  13c),3*10  one  might  expect  to  see  an  enhanced  phosphorescence  signal  in 
nitroxide-substituted  compounds  as  compared  with  their  diamagnetic  forms.  Observations  do 
not  bear  this  out  (Fig.  8);  triplet  (480  nm)  and  singlet  emission  bands  are  both  clearly  visible  for 
IA,  IB,  and  ID  at  77  K  (EPA  glass),  but  for  the  nitroxide  adducts  (IA,  ED)  the  triplet  emission 


is  diminished  to  a  greater  extent  than  that  of  the  singlet.  Thus  we  cannot  distinguish  between 
pathway  13c  and  a  rapid,  two-step  relaxation  via  13b. 

CoiifilUSiQDS 

For  all  adducts  examined  in  this  study,  absorption  and  fluorescence  emission  spectral 
shapes  and  energies  were  identical  to  those  of  the  parent  fluorophore.  There  was  no  evidence 
for  charge-transfer  absorption  or  charge-recombination  emission.  We  conclude  that  the 
emission  arises  from  the  same  locally-excited  singlet  state  in  the  presence  or  absence  of 
nitroxide  substituents. 

Intramolecular  quenching  of  excited  singlet  states  by  nitroxyl  radicals  is  very  efficient; 
rate  constants  are  in  the  range  of  lO^-lO1®  s'1.  At  distances  as  great  as  12  A,  the  quenching  can 
occurs  at  rates  >  108  s'*.  The  position  of  ring  substitution  and  the  type  of  linkage  appears  to 
significantly  influence  the  quenching  efficiency. 

Forster  energy  transfer  cannot  account  for  more  than  a  small  fraction  of  the  observed 
rates,  even  assuming  larger  values  for  x2.  The  rates  do  not  show  the  dependence  on  overlap 
integrals  predicted  by  theories  of  Forster  or  Dexter  energy  transfer.  There  is  no  clear  correlation 
of  kq  values  with  calculated  values  of  AGet.  suggesting  that  photo-induced  electron  transfer  is 
also  unimportant  Significant  quenching  rates  are  obtained  in  hematoporphyrin-nitroxide 
adducts  which  contain  a  fluorophore  whose  singlet  energy  is  below  that  of  the  lowest  excited 
state  of  the  nitroxide.  Since  the  driving  force  for  electron  transfer  is  also  small,  this  result 
suggests  that  the  nitroxide  moiety  facilitates  either  intersystem  crossing  to  the  triplet  or  internal 
conversion  to  the  ground  state.  Further  support  for  exchange-induced  relaxation  is  provided  by 
the  linear  correlation  between  knr  and  kq,  indicating  a  fixed  enhancement  of  a  pre-existing 
pathway.  Thus  we  conclude  that  the  quenching  arises  primarily  through  electron  exchange. 

This  conclusion  is  consistant  with  earlier  studies  of  diffiisional  quenching  of  aromatic 
hydrocarbons  by  nitroxides.lb 

This  work  provides  the  first  preliminary  evidence  for  long-range  quenching  by  electron 
exchange  that  proceeds  via  the  local  relaxation  of  the  singlet  and  not  energy  transfer.  That  this 


pathway  also  appears  to  outcompete  or  to  occur  in  preference  to  exergonic  electron  transfer  is 
somewhat  surprising,  given  current  thinking  on  the  relationship  between  electron  exchange  and 
electron  transfer  processes.  15b.30b.42  While  compounds  containing  fluorophores  with  different 
structures,  singlet  energies,  and  redox  potentials  could,  in  principle,  exhibit  a  number  of 
different  quenching  mechanisms,  data  from  this  and  previous  studies1  have  failed  to  provide 
unequivocal  evidence  for  the  operation  of  mechanisms  other  than  electron  exchange  in  the 
quenching  of  excited  singlet  states  by  nitroxides.  Further  work  employing  a  series  of  differing 
nitroxides  and  fluorophores  combined  within  a  rigid  structural  framework15  could  provide 
additional  insight  concerning  the  factors  controlling  the  possible  relaxation  pathways. 

Experimental 

Optical  Measurements 

Solutions  for  optical  measurements  in  all  solvents  were  prepared  from  concentrated 
(-10  mM)  stock  solutions  in  methanol.  All  solvents  were  HPLC  grade  (Aldrich)  used  without 
further  purification.  Absorption  spectra  were  measured  on  solutions  in  1  or  5  cm  quartz  cells 
using  a  Hewlett-Packard  8451 A  diode  array  spectrophotometer  (spectral  resolution,  2  nm). 

Quantum  yield  measurements  were  made  on  an  SLM-Aminco  SPF-500C 
spectrofluorometer  linked  to  an  IBM-XT  for  data  collection  and  processing.  Freshly  made 
solutions  had  optical  densities  of  <0.02  at  the  excitation  wavelength,  chosen  as  the  absorbance 
maximum  in  the  range  of  280-300  nm.  Excitation  and  emission  bandpasses  were  2  nm  in 
most  cases;  4  nm  bandpasses  were  used  to  measure  some  of  the  very  low  quantum  yields. 
Samples  were  thermostated  at  20*C.  Solutions  were  deoxygenated  through  bubbling  with  N2; 
emission  scans  were  repeated  until  bubbling  caused  no  further  increase  in  fluorescence 
intensity.  Any  signal  due  to  the  solvent  was  removed  by  subtraction.  Spectra  were  corrected 
for  instrumental  response  and  integrated  using  correction  factors  and  software  provided  by  the 
manufacturer.  (A  modified  correction  procedure  was  used  to  avoid  arbitrary  normalization  and 
conversion  steps  in  the  standard  routines). 
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Quantum  yields  were  calculated  with  respect  to  quinine  sulfate  in  0.1  N  H2SO4 

(OD  =  0.02)  and  corrected  for  refractive  index  differences  between  solvents,  using 

2 

,  F$  Aq$ns<}>qS 
(14)  <j>= - 

As  F qS  n0 

where  the  subscript  s  refers  to  the  sample,  qs  refers  to  quinine  sulfate,  $  is  quantum  yield,  A  is 
absorbance  at  the  excitation  wavelength,  F  is  the  integrated,  corrected  fluorescence,  ns  is  the 
refractive  index  of  the  sample  solvent,  and  no  is  the  refractive  index  of  water.  0qS  is  0.55. 

Fluorescence  lifetimes  were  measured  by  time  correlated  single  photon 
counting,  on  an  instrument  generously  made  available  through  the  Laser  Biomedical  Research 
Center  at  MIT.  Excitation  was  provided  by  a  cavity-dumped  tunable  dye  laser,  synchronously 
pumped  with  the  second  harmonic  of  a  mode-locked  Nd-YAG  (Coherent,  Antares  76-S).  For 
the  naphthalene  derivatives  the  dye  pulses  (5  ps  fwhm,  20  MHz)  were  frequency  doubled  with 
a  phase  matched  KTP  crystal  to  295  nm.  The  porphyrin  derivatives  were  excited  directly  with 
the  dye  laser,  tuned  to  580  nm.  Fluorescence  from  deoxygenated  samples  was  collected  at  90* 
(naphthalene  derivatives)  or  180'  (porphyrin  derivatives)  and  focused  into  a  Spex 
monochometer,  entrance  and  exit  slits  were  adjusted  to  obtain  a  counting  rate  of  4000-5000 
photons  per  second  impinging  on  the  PMT. 

A  Canberra  2143  time-amplitude  convener  (TAC)  received  start  pulses  directly  from 
the  cavity  dumper  and  stop  pulses  from  the  PMT.  The  TAC  signal  is  passed  to  a  multi-channel 
analyzer  and  the  data  are  saved  in  an  IBM-XT  for  later  analysis  with  a  Sun  microcomputer.  The 
instrumental  response  limit  for  this  system  is  about  70  ps. 


The  general  strategy  for  this  search  was  to  allow  the  functional  groups  of  each 
compound  to  freely  rotate  around  all  nonconju  gated  bonds  that  were  not  in  constrained  ring 
structures.  There  were  two  degrees  of  rotational  freedom  in  IA  and  IDA:  1)  the  bond  between 
the  carboxyl  carbon  and  oxygen  atoms  of  the  ester  linkage;  and  2)  the  carboxyl  oxygen  and 
carbon  Cl  of  the  nitroxyl  ring.  The  functional  groups  were  rotated  through  360*  around  each 


of  these  two  bonds  in  5’  increments.  The  structures  after  each  rotation  were  energy  minimized 
using  AMBER28  to  optimize  the  bond  distances,  bond  angles,  dihedral  angles  of  the  bonds, 
dipole  and  van  der  Waal's  interactions.  The  structures  were  allowed  to  minimize  to  a 
convergence  of  1.0  x  10*1  kcal/mol  difference  in  overall  energy,  and  0.10  kcal/mole  A  in  the 
gradient  of  the  overall  energy.  An  analysis  of  the  energies  associated  with  each  rotation  versus 
the  angle  of  rotation  around  each  degree  of  freedom  showed  that  the  local  energy  wells  were 
greater  than  +/-10*  in  width  relative  to  the  angles  of  the  minimized  conformations.  Thus  a 
smaller  rotational  increment  was  not  necessary  in  searching  for  the  global  energy  minima. 

The  addition  of  a  methylene  group  between  the  naphthalene  ring  and  the  carboxylate 
introduces  two  additional  degrees  of  rotational  freedom  to  HA.  To  simplify  this  system,  the 
structure  was  allowed  to  energy  minimize  in  terms  of  only  the  naphthalene  ring  and  the 
intervening  methylene  group.  This  fixes  the  methylene  hydrogens  in  orientations  that  are  not 
statically  hindered  by  the  aromatic  ring.  The  structure  was  then  treated  in  the  same  fashion  as 
IA  and  IIIA  through  rotations  around  the  two  carbon-oxygen  bonds  in  the  ester  linkage,  but  also 
including  the  additional  rotation  about  the  methylene-carboxylate  bond. 

Preparation  of  Nitroxide  Derivatives 

General 

All  starting  materials  and  reagents  were  obtained  from  Aldrich  and  used  as  received. 
High  purity  grade  Burdick  and  Jackson  solvents  were  used  following  distillation  over  CaH2  in 
all  preparations.  Chromatographic  separation  was  obtained  on  a  flash  column  using  30  micron 
silica  (Amicon  or  Davisil)  under  positive  N2  pressure,  solvents  are  indicated  where  appropriate. 
Melting  points,  which  are  uncorrected,  were  measured  on  a  Mel-temp  apparatus  (Laboratory 
Devices).  Infrared  (IR)  spectra  were  recorded  as  a  thin  film  on  a  Mattson  Sinus  100  infrared 
spectrometer  at  4  cm'*.  Mass  spectra  (MS)  were  obtained  by  thermal  desorption  from  a  copper 
probe,  in  either  an  electron  impact  (El)  or  chemical  ionization  (Cl)  mode  on  a  Finnigan  45 10 
spectrometer  interfaced  with  an  INCOS  2300  data  system.  Cl  mass  spectra  of  nitroxide 
derivatives  gave  an  (M  +  2)+  molecular  ion  suggesting  the  nitroxide  derivative  is  chemically 


reduced  prior  to  ionization  in  the  mass  spectrometer.  This  was  not  observed  in  the  spectra  of 
nitroxides  obtained  by  El  or  in  the  spectra  of  the  nitroxide-acetyl  derivatives  by  CL  Nuclear 
Magnetic  Resonance  (NMR)  spectra  of  the  acetyl  derivatives  were  obtained  on  either  a  Bruker 
AM-250  or  Varian  XL- 300  at  Massachusetts  Institute  of  Technology  and  measured  in  CDCI3. 
4-(2,2,6,6-tetramethylpiperidinyloxy,  free  radical)  1-naphthoate  (1A).  1- 
Naphthoyl  chloride  (950  mg)  was  dissolved  in  methylene  chloride  (10  ml)  and  pyridine  (1  ml). 
4-Hydroxyl-TEMPO  (860  mg  in  2.0  ml  of  methylene  chloride)  was  added  drop  wise  to  the 
stirred  solution.  After  20  h  this  solution  was  diluted  with  methylene  chloride,  washed  with 
water,  dried  (Na2S04)  and  evaporated  to  an  orange  oil  which  later  crystallized  on  standing. 
Column  chromatography  (BIO-RAD,  BIO-SIL  A  100-200  mesh)  eluting  with  methylene 
chloride  and  recrystallization  from  methanol  afforded  the  desired  product  as  orange  needles: 

962  mg;  mp  100-102*  C;  IR  1713  cm'1;  MS(EI)  326  (M+,  18),  155  (100),  MS(CI) 

328([M  +  2H]+,  25),  312  (M  -  14,35),  310(M  -  16,45),  156  (100). 
2-[4-(2,2,6,6-tetramethylpiperidinyloxy,  free  radical)]  6-hydrogen 
naphthylenedicarboxylate  (IVA)  and  2,6-Di[4-(2,2,6,6- 

tetramethylpiperidinyloxy,  free  radical)]  naphthylenedicarboxylate  (IVC).  2,6- 
Naphthalenedicarboxylic  acid  (3.5  g)  was  suspended  in  methylene  chloride  (100  ml)  and 
pyridine  (3.0  ml).  Thionyl  chloride  (3.0  ml)  was  then  added  dropwise  and  the  resultant 
solution  was  stirred  until  no  suspension  remained  and  the  solution  was  clear.  The  intermediate 
diacid  chloride  was  evaporated  to  dryness  in  an  effort  to  remove  extraneous  thionyl  chloride 
then  dissolved  in  methylene  chloride  (150  ml)  and  pyridine  (3.0  ml).  4-Hydroxy-TEMPO  (1  g) 
in  a  small  amount  of  methylene  chloride  was  added  and  this  mixture  stirred  overnight.  Water 
(3.0  ml)  was  added  and  the  precipitated  diacid  was  filtered  off  and  washed  with  excess 
methylene  chloride.  The  combined  filtrate  was  washed  with  water,  dried  with  anhydrous 
sodium  sulfate  and  evaporated  to  a  dark  orange  oil.  The  desired  products  were  separated  and 
purified  by  flash  chromatography  (eluting  with  3%  methanol/methylene  chloride)  and 
recrystallized  from  methanol  giving  compound  IV  A  (90  mg)  and  compound  IV  D  (220  mg). 
Compound  IVA:  mp  >  200*  C(dec);  IR  3200-2200  cnr1  (b),  1715, 1622  cm*1;  MS(EI) 


370(M+),  199(80),  124(100).  Compound  IVD:  mp  >  225*  C(dec);  ER  1712  cm'1;  MS(EI) 
524(M+,  30),  370(20),  124(100). 

4-(2,2,6,6-tetramethyipiperidinyloxy,  free  radical)  2-naphthoate  (IIIA).  2- 
Naphthoyl  chloride  (380  mg)  and  4-Hydroxy-TEMPO  (345  mg)  were  dissolved  in  cold 
pyridine  and  stirred  overnight  (20  hr).  This  mixture  was  diluted  with  methylene  chloride, 
washed  with  water  (1  x  25  ml),  1  HC1  (3  x  25  ml),  dried  over  an  hydrous  sodium  sulfate, 
and  evaporated  to  dryness.  The  desired  product  was  purified  by  flash  chromatography 
(methylene  chloride),  obtained  as  a  solid  and  recrystallized  from  methylene  chloride/n-hexane: 
235  mg;  mp  122-123*  C;  IR  1713  cm'1:  MS(CI)  328  [M  +  2]+  75),  312(M  -  14,70), 

310(M  -  16,50),  156(100). 

4-(2,2,6,6-tetramethylpiperidinyloxy,  free  radical)  (1-naphthyl)  acetate  (IIA). 
1-Naphthylacetic  acid  (500  mg)  was  combined  with  dicyclohexylcarbodiimide  (DCC,  610  mg), 
dimethylaminopyridine  (DMAP,  32  mg)  and  4-Hydroxy-TEMPO  (462  mg)  and  dissolved  in 
methylene  chloride  and  stirred  overnight  under  an  atmosphere  of  N2^.  The  solid  dicyclohexyl 
urea  was  filtered  off  and  further  rinsed  with  methylene  chloride.  The  combined  filtrate  was 
washed  with  IN  HC1  (2  x  75  ml),  water  (3  x  75  ml),  dried  (Na2S04)  and  evaporated.  The 
product  was  purified  by  flash  chromatography  (methylene  chloride)  and  evaporated  to  an  orange 
oil  which  did  not  crystallize:  IR  1733  cm"1  MS(CI)  342([M  +  2]+,  25),  326(M  -  14,40), 
324(M  -  16,38),  156(100). 

4-[17p-hydroxy-4,4'-dimethyIspiro  (5a-androstane-3,2'-oxazolidin)-3'yloxy, 
free  radical]  l-naphthoate(ID).  Compound  ID  was  prepared  by  the  DCC  coupling 
described  above  using  1-Naphthoic  acid  and  3-DOXYL-17(3-hydroxy-5  a-androstane 
(Aldrich).  The  desired  product  was  purified  by  TLC  on  a  1  mm  thick  plate  (20  x  20  cm, 
Whatman)  eluted  with  methylene  chloride  and  obtained  as  a  yellowish  solid:  mp  153  -  154*  C; 
IR  1714  cm’1;  MS(EI)  530(M+,  1),  444(DOXYL  fragment  loss,  2),  155(100). 

Preparation  of  Nitroxide-Acetvl  Derivatives 

4-(2,2,6,6-tetramethylpiperidinylacetoxy)  1-naththoate  (IB).  Compound  IA  (50 
mg)  dissolved  in  tetrahydrofuran  (3  ml)  in  the  presence  of  5%  palladium/carbon  (5  mg)  was 


stirred  under  an  atmosphere  of  hydrogen  for  3  h.  After  this  time  the  solution  was  filtered 
directly  (with  minimal  exposure  to  air)  into  a  solution  of  acetic  anhydride/pyridine  (1:10)  which 
had  been  previously  deoxygenated  with  nitrogen.  This  solution  was  stirred  for  20  h  under  N2 
then  diluted  with  methylene  chloride,  washed  with  water,  dried,  and  evaporated  to  a  yellowish 
oil.  After  chromatography  on  silica  (eluting  with  methylene  chloride)  fractions  containing 
product  crystallized  after  solvents  were  removed.  The  final  product  was  recrystallized  from 
methanol  to  remove  a  small  amount  of  compound  IA  and  gave  a  snowy  white  solid:  20  mg;  mp 
86-87*  C;  IR  1712,  1770  cm*1;  MS(EI)  369(M+),  327(M+  -  42),  155(90),  140(100); 
NMR(250  MHz)  8  8.80(d),  8.14(d),  8.02(d),  7.88(d)  7.55(m)  (7H,2,3,4,5,6,7,8  -H),  5.45 
(m,  11  -H),  2.13(s,OCOMe),  2.20  -  2.00(m,  4H,  12  -CH2),  1.35,  1.17(s,  12H,  13a  -Me). 
2-[4-(2,2,6,6-tetramethylpiperidinvIacetoxy)]  6-hydrogen 
naphthylenedicarboxylate  (IVB).  2,6-Di[4-(2,2,6,6- 

tetramethylpiperidinylacetoxy)]  naphthylenedicarboxylate  (IVE),  4-(2, 2,6,6- 
tetramethylpiperidinylacetoxy)  2-naphthoate  (IIIB),  4-(2, 2,6,6- 
tetramethylpiperidinylacetoxy)  (l-naphthyl)acetate  (IIB).  These  compounds  were 
prepared  in  a  fashion  analogous  to  compound  IB.  However,  the  hydroxylamine  produced  in 
the  catalytic  reduction  of  compound  IVA  precipitated  upon  its  formation  and  excess 
tetrahydrofuran  and  pyridine  was  added  in  order  to  filter  off  the  catalyst  Compound  IVB  was 
not  purified  chromatographically  but  was  obtained  as  a  white  solid  from  crystallization  in 
methanol.  Compound  IVB:  IR  1769, 1753,  1715  cm*1;  MS(EI)  413(M+),  371  (M+  -  42), 
140(100);  NMR(250  MHz)  8  8.74,  8.61  (s,l,5  -H),  8.22(d),  8.18(m)(4H  3,4,7,8  -H), 
5.45(m,  12  -H),  2.16(s,  OCOMe),  2.11(m,  4H,  13  -CH2),  1.34,  1.18(s,  12H,  14 -Me). 
Compound  IVE:  IR  1766,  1703  cm*1;  MS(EI)  610(M+,5),  568(M+  -42,30),  124(100); 
NMR(250  MHz)  8  8.58(s,  1,5 -H),  8.04(m,3,4,7,8  -H),  5.43(m,12 -H),  2.13(s,  OCOMe), 
2.03(m,4H,  13 -CH2),  1.33,  1.16(s,  12H,  14 -Me). 

Compound  mB:  mp  115-116*  C  IR  1732  cm*1,  MS(CI)  398  (M  +  29,5),  370([M  +  1]+, 
100),  310(M  -  CH3C02H.95),  198(30),  NMR(250  MHz)  8  8.54(s,  1H,  1-H),  8.00(d), 


7.94(d),  7.86(d),  7.55(m)  (6H,  3, 4,5, 6,7,8  -H),  5.39(m,  11  -H),  2.1  l(s,OCOMe),  2.15  - 
1.90  (m,  4H,  12-CH2),  1.31,  1.14(s,  12H,  13a-Me). 

Compound  nB:  IR  1767,  1732  cm'1,  MS(CI)  412(M  +  29,5),  384([M  +  1]+,  100),  324 
(M  -  CH3CO2H,  100),  198(45),  NMR(250  MHz)  5  7.96(d),  7.84(d),  7.77(d),  7.49(m), 
7.38(m)(7H,  2,  3,  4,  5,  6,  7,  8  -H),  5.06  (m,  11  -H),  4.02(s,  2H,  la-CH2),  2.06(s, 

OCOMe),  1.90-  1.67(m,  4H,  12 -CH2),  1.17,  1.04(s,  12H,  13a-Me). 

Preparation  of  Methvl  Esters 

Methyl  1-naphthoate  (IC)  and  Methyl  2-naphthoate  (IIIC).  1-Naphthoyl  chloride 
(500  mg)  1  or  2-Naphthoyl  chloride  was  dissolved  in  cold  pyridine  (3  ml),  methanol  (2  ml)  was 
added,  and  the  resultant  mixture  was  stirred  in  an  ice  bath  for  2  hr  and  at  room  temperature  for 
an  additional  4  hr  after  which  the  reaction  mixture  was  diluted  with  methylene  chloride,  washed 
with  water  (1  x  25  ml),  1M  HC1  (3  x  25  ml),  water  (2  x  25  ml) ,  dried  (Na2SC>4),  and 
evaporated  to  dryness.  The  product  was  purified  by  flash  chromatography  (methylene 
chloride).  Compound  IC:  obtained  as  a  colorless  oil;  IR  1711  cm'1.  Compound  IIIC: 
obtained  as  a  solid  and  recrystallized  from  methylene  chloride/n-hexane  mp  72  -  73*  C,  IR 
1711  cm'1,  MS(CI)215([M  +  29]+,  15),  187([M+  1]+,  100),  155(10). 

2-Methyl  6-hydrogen  naphthylenedicarboxylate  (IVC),  2,6-dimethyl 
naphthodylenedicarboxylate  (IVF)  and  Methyl  l-(naphthyl)acetate  (IIC). 
Compounds  IVC  and  F  were  prepared  by  the  procedure  used  in  the  synthesis  of  compounds 
IVA  and  D  substituting  methanol  for  4-Hydroxy-TEMPO.  Compound  IVC:  mp  >  200’  C 
(dec.),  IR  1721,  1682  cm'1,  MS(CI)259(M  +  29,15),  231([M  +  1]+,  100),NMR  (250  MHz, 
CDCL3  and  d5.pyridine)  5  8.69,  8.6I(s,l,4  -H),  8.20(d  of  d),  8.08(d  of  d),  7.98(d)  (4H,  3, 

4,  7,  8  -H),  3.96  (s,  3H,  2a-OMe).  Compound  IVF:  mp  188*  C  (lit.  186’  C),  IR  1707  cm'1. 
Compound  IIC  was  prepared  by  the  procedure  used  for  the  synthesis  of  compound  IIB 
substituting  methanol  for  the  nitroxide:  obtained  as  a  colorless  oil;  IR  1736  cm'1;  MS(CI) 
229([M  +  29]+,  8),  201([M  +  1]+,  100),  169([M  +  1]+  -OMe,  65). 

Preparation  of  Hematoporphvrin  DC  carboxv-PROXYL  Derivatives 


2(4)-Vinylhematoporphyrin  IX  DME  was  obtained  from  commercial  hematoporphyrin  EX 
(Aldrich)  by  the  H2SC>4/methanol  procedure  reported  by  Smith.38  The  product  was  purified 
by  flash  column  on  silica  (3  -  5%  methanol/methylene  chloride)  and  obtained  as  a  solid  from 
methylene  chloride/n-hexane.  At  this  point,  no  effort  was  made  to  resolve  the  isomeric  mixture 
of  2(4)-Vinylhematoporphyrin  IX  DME. 

Vinvlhematoporphvrin  IX  DME  IR  3442  cm‘^(bd,  O  -  H),  3319  (m,  N-H),  1734(s,  ester 
C=0);  MS(CI)637([M  +  29]+,  8),  609([M  +  1]+,  70),591([M  +  1]+  -H2O,  100);  NMR 
(300  MHz)  8  10.18,  10.06,  10.01,  9.93,  9.91,  9.88  (4H,a,b,g  and  d  meso  -H),  8.20  (m,  1H, 
2a  and  4a  -H  of  vinyl),  6.35,  6. 17  (m,  3H,  2a  and  4a  -H  of  hydroxyl  ethyl  and  2b  and  4b 
-CH2),  4.31  (4H,  6a  and  7a-CH2),  3.65,  3.61,  3.55,  3.53,  3.50,  3.43,  3.42  (overlapping  s's, 
18H,  1 ,3,5,8  -Me  and  6,7  -OMe),  3.23  (t,  4H,  6b  and  7b  -CH2),  2.12  (d,  3H,  2a  and  4a-Me). 
2(4)-l[-(hydroxy)ethyl]-4(2)-ethyl-6,7-bis[2-(methoxycarbonyI)ethyl]-l,3,5,8- 
tetramethyl  porphyrin  (VC) 

2(4)-Vinylhematoporphyrin  IX  DME  dissolved  in  *- "-ahydrofuran  in  the  presence  of  5% 
palladium/carbon  (5  mg)  was  stirred  under  an  atmosphere  of  hydrogen  overnight  After 
filtering  off  the  catalyst  the  desired  product  was  purified  on  a  flash  column  (3% 
methanol/methylene  chloride)  and  obtained  as  a  solid:  IR  3454  cm‘l  (bd,  O  -  H),  3318  (m, 

N  -  H),  1735  (s,  ester  C=0);  MS(CI)  639  «M  +  29]+,  10),  61 1([M  +  1]+,  100), 

593([M  +  1]+  -H2O,  95);  NMR(300  MHz)  5  10.34,  10.03,  10.01,  10.00,  9.99, 
9.97(4H,a,b,g  and  d  meso  -H),  6.32(q,  1H,  2a  and  4a  -H  of  hydroxy  ethyl),  4.34  (m,  4H,  6a 
and  7a-CH2),  4.05  (m,2H,2a  and  4a  -CH2),  3.66,  3.65,  3.64,  3.57,  3.56,  3.55  (overlapping 
s's,  18H,  1,3,5,8-Me  and  6,7-OMe),  3.25(t,4H,6b  and  7b-CH2),  2.14(d,3H,2a  and  4a-Me  of 
hydroxy  ethyl),  1.84(t,3H,2a  and  4a- Me  of  ethyl). 
2(4)-(l-[3-(2,2,5,5-tetrainethylpyrrolidinyIoxy,  free 
radical)carbonyloxy])ethyI-4(2)*ethyl-6,7*bis[2-(methoxycarbonyl)ethyl]< 

1,3,5, 8-tetramethylporphyrin  (VA,B).  The  nitroxide  derivative  was  prepared  by  the 
DCC  coupling  reaction  with  3-carboxy-PROXYL  and  the  hydroxy  porphyrin.  A  and  B  isomers 
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were  separated  by  TLC  (3%  methanol/methylene  chloride)  but  were  not  uniquely  identified  as  to 
substituant  position  (2  or  4). 

2(4)-carboxv-PRQXYL-4(2)-ethvlhematoporphvrin  IX  DME  (VA):  IR  3316  crn'l;  (m,N  -H), 

1735(s, ester  c=0);  MS(CI)807([M  +  29]+,  1),  779([M  +  l]+,3),  764([M  +  1]+  -15,8), 

593([M+ 1J+-H20,  100). 

Supplementary  Material  Available:  A  listing  of  atomic  coordinates  for  low  energy 
conformations  of  IA,  13A  and  IIIA  (8  pages).  Ordering  information  is  given  on  any  current 
masthead  page. 
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Figure  Captions 

Figure  1.  Naphthalene  derivatives  examined  in  this  study. 

Figure  2.  Fluorescence  excitation  and  emission  spectra  of  IYD  in  hexane. 

Figure  3.  Fluorescence  decays  of  nitroxide  adducts  (b)  as  compared  to  diamagnetic  analogues 
(a):  (1)  ID  and  IB  in  acetonitrile,  (2)  ID  and  IB  in  hexane,  (3)  IIA  and  IIB  in  acetonitrile, 

(4)  I  IIA  and  DIB  in  acetonitrile.  Excitation  wavelength  was  295  nm.  Emission  was  monitored 
at  the  peak  maximum  for  each  compound  (see  Table  1). 

Figure  4.  Absorption  spectra  of  TEMPOL  in  a  series  of  solvents  showing  the  solvent-dependent 
shifts  in  the  absorption  maximum  and  extinction  coefficients  (M*1  cm'1).  Fluorescence  emission 
spectrum  of  IB  in  methanol  is  included  in  order  to  illustrate  the  degree  of  spectral  overlap 
between  naphthalene  emission  and  nitroxide  absorption. 

Figure  5.  Intramolecular  quenching  constants  vs  spectral  overlap  integrals  for  IA(B),  IIA(A), 
and  IIIA(#);  (a)  Forster  spectral  overlap  integral,  Jf  (M  'cm3).  The  rate  constants  have  been 
normalized,  k'=kq  Tdn4/<j><j,  to  account  for  solvent  dependent  variations  in  lifetime,  quantum 
yield  and  refractive  index,  (b)  Dexter  spectral  overlap  Jd  (cm). 

Figure  6.  Hematoporphyrin  derivatives  examined  in  this  study.  The  two  isomers  (V  A  and  B) 
with  the  nitroxide  at  either  the  2  or  4  position  were  separated,  but  were  not  uniquely  identified 

Figure  7.  Intramolecular  quenching  constants  (kq)  of  nitroxide  adducts  IA(B),  IIA(A), 
IHA(#),  IVA(A),  IVD(X),  and  VA(*)  vs  non-radiative  rate  constants  (knr)  of  the 
corresponding  diamagnetic  O-acetyl  and  methyl  ester  derivatives,  kq  values  for  IVD  are  divided 
by  two  to  account  for  the  presence  of  two  nitroxide  groups  in  this  compound.  The  line 
represents  a  linear  regression  (slope  =  0.96,  r2  =  0.76)  of  all  data  excluding  IVA  in  methanol 
(circled)  and  VA.  knr  for  VA  in  methanol  was  estimated  from  data  of  Andreoni,  et.al.43 

Figure  8.  Fluorescence  (1)  and  phosphorescence  (2)  emission  spectra  of  IA,  IB,  and  ID. 
Excitation  wavelength,  300  nm;  bandpass,  2  nm;  concentration,  15  (iM  in  EPA  glass,  77  K. 
Phosphorescence  intensities  are  expanded  by  a  factor  of  5. 
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Table  1.  Absorption  and  fluorescence  maxima  of  MV1 

absorbance  in  methanol  fluorescence  emission  maxima  (nm) 


compounds 

abs  max  (nm) 

£  (M’  cm'1) 

hexane 

acetonitrile 

methanol 

water 

Series  1 

222 

4.2  x  104 

332 

361 

30? 

390 

300 

7.1  xIO3 

343 

358 

Series  II2 

228 

- 

328 

326 

326 

325 

274 

340 

335 

337 

335 

284 

318 

353 

294 

352 

371 

Series  III 

242 

1.9  x  104  (8.9  x  104) 

336 

344 

364 

376 

276 

7.9  xIO3 

352 

359 

351 

284 

9.6  x  103  (1.3X104) 

367 

375 

296 

6.8  xIO3  (9.7  xIO3) 

388 

324 

1.6 xIO3  (2.2 xIO3) 

338 

1.9 xIO3  (2.7 xIO3) 

Series  IV  A-C 

248 

7.0  xIO4 

355 

359 

359 

360 

288 

1.2  xIO4 

373 

377 

374 

381 

298 

1.1  xIO4 

393 

397 

396 

334 

1.8  xIO3 

416 

426 

348 

2.0  xIO3 

IV  Df 

246 

8.3 xIO4  (8.4 x  104) 

354 

360 

363 

286 

1.4  xIO4 

372 

378 

380 

296 

1.4  x  104  (1.3  x  104) 

391 

400 

338 

2.1  xIO3  (2.3 xIO3) 

425 

352 

2.4 xIO3  (2.8 xIO3) 

Absorption  spectra  and  extinction  coefficients  measured  in  methanol.  Values  in  parentheses 
are  extinction  coeffients  for  methyl  esters  (III  C  and  IV  F)  for  which  the  coefficient  differs  slightly 
from  that  of  the  nitroxide  derivative.  The  nitroxide  contributions  to  the  absorption  of  l-IVA  and 
UVD  are  described  in  the  text  and  in  fig.  4. 

2Series  II  compounds  were  isolated  as  oils;  extinction  coeficients  were  not  measured. 


Table  2.  Fluorescence  lifetimes  and  quantum  yields  of  l-IV*. 
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Table  3.  Rate  constants  for  Intramolecular  quenching  in  the  naphthalene-nitroxide  adducts1 
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Table  4.  Total  energies,  nitrogen  atom-naphthalene  ring  distances  and 
k2  for  lowest  energy  conformations  of  IA,  ID,  IIA,  and  IIIA. 


compound 

energy1 

distance2 3 4 5 6,  R 

K23 

kcal/mol 

• 

A 

IA 

2225 

824 

0.18 

22.30 

826 

0.11 

22.39 

8.25 

012 

average4 

22.31 

8.25 

0.14 

—  ID5 

12.2 

IIA 

14.64 

6.55 

0.01 

15.01 

7.81 

0.03 

15.02 

7.41 

0.16 

15.18 

7.94 

0.07 

15.42 

8.32 

0.21 

15.47 

8.66 

2.01 

15.54 

8.83 

0.01 

10SQ 

152 

on 

average6 

15.12 

7.59 

0.26 

IIIA 

16.67 

9.04 

0.21 

16.68 

8.94 

0.18 

i&z 4 

9.04 

0.33 

average4 

16.70 

9.01 

0.24 

1  Energy  ror  each  conformation  is  the  sum  of  contributions  for  bond  lengths  and  angles, 
dihedral  bond  angles,  and  dipole  and  van  der  Waal's  interactions. 

2  ^stances  were  measured  from  the  nitrogen  atom  to  the  center  of  the  C9-C10  bond  of 
the  naphthalene  ring. 

3  k2  was  calculated29  from  the  emission  (naphthalene)  and  absorption  (nitroxide) 
dipoles.  Emission  dipoles  for  1-  and  2-substituted  naphthaldehydes  C1!^  transition)  were 
obtained  from  reference  40.  The  absorption  dipole  for  the  nitroxide  group  was  taken  to 
be  perpendicular  to  the  NO  bond  and  lying  in  the  CNC  plane,  based  on  the  local  Cs 

symmetry  of  the  nitroxide  group  and  the.analogy  to  the  n  —>n  *  transition  in 
carbonyls.15*1-22-23'41  0  <  k2  <4,  allowed  range. 

4  Numerical  averages  of  the  energies,  distances  and  k2  for  the  conformations  of  each 
compound. 

5  Distance  for  ID  was  estimated  from  molecular  models  using  the  same  relative 
orientation  of  the  carbonyl  group  as  that  calculated  for  the  lowest  energy  conformation 
of  IA. 

6  Averages  for  IIA  were  weighted  according  to  the  population  in  the  lowest  energy 
conformation  (E=14.64  kcal/mol)  compared  with  the  other  seven  states  (Eaverage=15-35 
kcal/mol).  The  ratio  of  population  in  states  2-8  relative  to  state  1  is  estimated  by:  N@- 
8)/N(i)=g(2-8)/0  (1)X  exp(-AE/RT>  7  exp(-(  15.35-14.64)/RT>  2.11. 
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Table  5.  Calculated  FSrster  energy  transfer  rates1. 


kpr x  lO9  s_1 


compound 

hexane 

acetonitrile 

methanol 

water 

A 

r  6.016  “ 

6.014 

6.150 

IIA 

0.006 

0.004 

0.008 

0.014 

IIIA 

0.010 

0.009 

0.024 

0.060 

N/A 

0.010 

0.026 

0.025 

0.025 

1  kpT  calculated  from  equation  7a  using  the  average  values  of  R  and  k2  (Table  4)  and 
<t>d.  Td  values  from  Table  2;  Jf  values  generated  from  equation  6a. 
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Table  6  Singlet  energies  (E^1  and  calculated  thermodynamic  driving  forces  (a  G£T)  for 
electron  transfer.2 


hexane3-4  acetonitrile2  water6 


compound 

^00 

a  Get 

a  Get 

Em 

aGet 

IA 

3.82 

-0.25 

3.43 

-0.90 

3.18 

-1.03 

ID^ 

3.82 

+0.03 

3.43 

-0.90 

IIA 

4.00 

+0.07 

4.00 

-0.87 

4.00 

-1.85 

IIIA 

3.72 

-0.07 

3.65 

-1.12 

3.30 

-1.14 

IVA 

3.54 

-0.29 

3.51 

-1.33 

3.44 

-1.28 

IVD 

3.54 

-0.29 

3.52 

-1.34 

. 

1  Singlet  energies  (Eoo)  were  obtained  experimentally  from  the  intersection  points  of 
fluorescence  excitation  and  emission  spectra. 

2  Calculated  from  equation  8  in  text;  units  are  eV. 

3  Favored  charge-separated  product  is  (Nap*-NO+).  Oxidation  potential  of  tempol 
(+0.63  V  vs.  SCE)  in  acetonitrile  from  reference  2;  reduction  potentials  for  1- 35 , 2- 351x0 
and  2.6- 350  substituted  naphthoates  (- 1 .92  V,  - 1 .95  V,  and  - 1 .55  V.  respectively ),  and 
l-methyl-naphthalene35b  (-2.5  V)  in  dimethylformamide,  acetonitrile  or  75% 
dioxane/water  were  employed  to  calculate  A  Go.  for  IA,  IIIA,  IVA,  and  IIA  respectively. 

4  Values  in  hexane  adjusted  for  solvation  (see  text)  using  eqn  10  and  the  approximation 
fd=fd=3.8  A . 

5  In  water  the  favored  charge-separated  product  is  (Nap+-NO*).  The  reduction  of  tempo 
gives  an  irreversible  wave  in  cyclic  voitametry  (-0.62  V)-32,33  The  oxidation  potentials  of 
l-35b  and  2- 354:3  methyl-naphthalenes  (+1.53  V  and  +1.55  V,  respectively)  in  acetic  acid 
solution  were  used  as  the  best  available  estimates  for  potentials  of  compounds 
substituted  at  the  1  and  2  positions. 

5  Redox  potentials  for  DOXYl  were  unavailable;  the  values  for  tempol 
were  assumed. 
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Synopsis 


Polarized  electronic  absorption  spectra  of  the  (100)  face  of  single  crystals  of  the  Z-form 
double  helical  duplex  of  d(m5CGUAm5CG)  have  been  obtained  from  Kramers-Kronig  analysis 
of  reflection  data.  The  c  crystallographic  axis  is  parallel  to  the  helix  axis  and  shows  but  weak 
absorption.  The  b  axis  is  perpendicular  to  the  helix  axis  and  shows  a  structureless  absorption 
band  centered  at  270  nm  with  an  oscillator  strength  of  0.26.  Calculations  of  the  crystal  spectra 
utilizing  available  transition  moment  data  for  the  individual  chromophores  are  carried  through 
using  the  oriented  gas  model  (no  interbase  interactions)  and,  again,  employing  all  base-base 
interactions  (point  dipole)  in  the  duplex.  The  calculated  hypochromism  of  the  270  nm  band  is 
much  less  than  the  experimental  value  obtained  from  the  crystal  data.  The  crystal  spectra  appear 
to  be  representative  of  Z-form  double  helices  of  essentially  infinite  length  and  not  of  a  collection 
of  twelve  base  duplexes.  No  evidence  for  nrc*  transitions  polarized  parallel  to  the  helix  axis  is 
found. 


INTRODUCTION 

At  the  simplest  level  of  interpretation,  the  electronic  absorption  spectrum  of  a  polynucleo¬ 
tide  is  just  a  summation  of  the  absorption  spectra  of  the  individual  chromophores  that  are 
present.  Refinements  bring  changes  owing  to  differences  in  the  environment  of  the  chromo¬ 
phores  and,  as  a  separate  issue,  chromophore-chromophore  interactions  (both  base  stacking  and 
pairing).  These  effects  lead  to  conformationally  dependent  changes  in  the  optical  properties  of 
solutions  of  polynucleotides  of  all  strandedness.  The  hypochromism  of  native  DNA  and  other 
stacked  polynucleotides  relative  to  unstacked  or  random  coil  arrangements  is  routinely  used  to 
assay  the  amount  of  order  present.1  More  recently  the  observation  that  Z-form  DNA  exhibits  less 
integrated  absorption  in  the  260  nm  region  than  even  the  B-form  has  been  used  to  study  various 
aspects  of  Z-B  transformations.2  Although  there  is  overall  hypochromism  exhibited  by  the  260 
nm  absorption  band  in  the  B  to  Z  transition  there  is  in  fact  pronounced  hyperchromicity  in  the 
long  wavelength  edge  (-295  nm). 
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The  absorption  curves  of  isotropic  solutions  are  of  course  a  spherical  average  of  all  orienta¬ 
tions  of  the  statistical  distribution  of  conformations  present,  and  since  the  broad,  diffuse  bands 
resulting  from  the  several  transitions  of  the  individual  chromophores  are  substantially  over¬ 
lapped,  little  detailed  electronic  spectral  information  is  available  from  solution  absorption  spec¬ 
tra.  Considerable  effort  has  been  applied  to  the  measurement  and  interpretation  of  CD  spectra  in 
order  to  elucidate  the  structure  and  conformation  of  polynucleotides.  Although  CD  spectra  are 
very  sensitive  to  conformation  and  contain  the  extra  dimension  of  the  sign  of  the  CD  bands,  it  is 
still  not  a  straightforward  matter  to  use  CD  data  to  gain  specific  information  regarding  the  effects 
of  base  pairing  and  stacking.  Linear  dichroism  (both  flow-  and  electric  field-induced)  studies  are 
useful  in  identifying  the  strengtii  of  absorption  parallel  and  transverse  to  the  helical  axis,  but 
necessary  estimates  concerning  the  extent  of  order  tend  to  limit  the  results  quantitatively,  and,  at 
best  there  is  still  a  random  distribution  about  the  orientation  axis. 

Recent  advances  in  the  crystallography  of  oligonucleotide  single  crystals  affords  us  the 
opportunity  to  examine  anisotropies  that  may  exist  in  the  absorption  tensor  of  these  polymers  by 
studying  the  polarized  electronic  spectra  of  such  crystals.  It  may  then  be  possible  to  relate 
observed  anisotropies  to  specific  conformations,  base  sequences  and  interactions  and  to  provide 
a  useful  test  for  the  theories  dealing  with  the  optical  properties  of  polynucleotides.  We  hope  that 
the  present  work  on  perfectly  ordered  systems  (viz.  crystals)  will  complement  the  other  experi¬ 
mental  approaches  mentioned  above. 

In  this  paper  we  report  polarized  UV  spectra  taken  from  the  be  face  of  a  single  crystal  of 
the  Z-form  d(m5CGUAm5CG)-d(m5CGUAm5CG)  duplex.  Since  the  c  crystallographic  axis 
coincides  with  the  helix  axis,  we  obtain  a  spectrum  for  radiation  polarized  parallel  to  the  helix 
axis  and  another  spectrum  for  radiation  polarized  perpendicular  to  the  helix  axis.  The  latter 
spectrum  (6-axis  polarized)  represents  the  transitions  of  the  duplex  as  they  project  onto  the  b 
axis  and  may  be  different  from  the  corresponding  projections  onto  the  a  crystallographic  axis.  It 
is  just  such  a  difference  that  we  would  ultimately  hope  to  exploit  when  the  a  axis  spectrum 
becomes  available.  For  comparison  we  report  spectra  for  all  thre-  crystal  axes  that  are  calculated 
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from  the  known  spectral  parameters  of  the  purine  and  pyrimidine  chromophores  using  an 
oriented  gas  approach  (no  interactions  or  mixing)  and  a  point  dipole  treatment  including  all 
base-base  interactions  of  the  twelve  member  duplex. 

Crystal  Structure 

Single  crystals  of  the  self-complementary  d(m5CGUAm5CG)  oligomer  are  of  orthorhombic 
space  group  P212121  and  have  unit  cell  dimensions  a  =  17.52  A,  b  =  30.44  A  and  c  =  44.82  A.3 
There  are  four  duplexes  of  the  Z-conformadon  per  unit  cell.  Duplexes  pack  end-to-end  not  only 
as  pairs  in  a  given  unit  cell  but  in  adjacent  cells  in  such  a  way  as  to  create  uniformly  continuous, 
parallel  strands  of  Z-DNA  with  the  helix  axis  parallel  to  the  crystallographic  c  axis.  Each  duplex 
makes  one-half  turn  of  the  helix.  Figure  1  shows  a  projection  of  six  (stick  model)  duplexes  onto 
the  be  face.  The  duplex  at  site  2  (x  +1/2,  -y,  z  +  1/2)  continues  the  half  helix  of  site  1  (x,y,z)  to 
form  one  full  turn  of  the  helix  in  one  unit  cell.  Projections  of  the  remaining  two  duplexes  of  the 
unit  cell  at  site  3  (x,  y  +  1/2,  z  +  1/2)  and  site  4  (x  +  1/2,  y  +  1/2,  z)  are  obtained  from  sites  1  and  2 
by  a  translation  and  then  a  180°  rotation  around  an  axis  parallel  to  b  (two- fold  screw).  Thus,  the 
transition  moment  of  a  particular  band  of  any  given  residue  at  all  four  unit  cell  sites  project 
equivalently  (except  for  phase)  onto  the  three  crystallographic  axes.  Since  observed  absorption  is 
proportional  to  the  square  of  these  projections,  all  four  duplex  sites  therefore  are  spectroscopi¬ 
cally  equivalent  when  viewed  along  any  unit  ceil  axis.  The  normals  to  the  planes  of  the  bases 
are  pitched  at  an  average  of  8°  from  the  c  axis.  The  direction  cosines  of  these  normals  and  the 
C4-C5  reference  axis  of  the  purines  and  the  N1-C4  reference  axis  of  the  pyrimidines  for  the 
twelve  sites  are  given  in  Table  L  The  attitudes  of  all  the  bases  in  the  structure  are  completely 
specified  by  the  values  given  in  Table  I. 

Another  presentation  of  the  crystal  structure  is  shown  in  Fig.  2,  where  the  atoms  are  drawn 
with  their  van  der  Waal’s  surfaces.  The  heavy  solid  lines  connect  the  phosphorus  atoms  of  the 
phosphate  backbone.  The  packing  of  the  duplexes  into  continuous  helical  strands  with  no  "gaps" 
is  apparent  in  this  figure.  Crystals  were  grown  by  a  vapor  diffusion  method  that  will  be 
described  elsewhere.3  Relatively  large  crystal  faces  (-200  |i)  are  needed  for  the  spectral 


measurements,  and  crystals  of  this  size  were  occasionally  obtained.  The  general  morphology  is 
given  in  Fig.  3.  The  crystals  are  pseudohexagonal  and  elongated  along  c.  The  principal  faces  are 
(100),  (001)  and  (110). 

Spectra 

Owing  to  the  orthorhombic  crystal  symmetry,  the  be  crystal  face  will  show  two  "indepen¬ 
dent"  absorption  spectra.  These  spectra  correspond  to  the  incident  radiation  polarized  parallel  to 
the  b  axis  and  parallel  to  the  c  axis  and  will  sample  transitions  to  B2  and  Bj  crystal  factor  group 
states  respectively.  Radiation  polarized  at  intermediate  angles  will  simply  be  a  cosine  squared 
weighted  average  of  the  b  axis  and  c  axis  spectra.  Based  on  general  considerations  of  the  nature 
of  the  electronic  spectra  of  the  purine  and  pyrimidine  chromophores  one  expects  relatively 
strong  absorption  polarized  perpendicular  to  the  helix  axis  (7t-*rc\  in-plane  polarized)  and  about 
10  to  100  times  less  absorption  polarized  parallel  to  the  helix  axis  (vide  infra).  The  molar  con¬ 
centration  of  base  chromophore  is  3.30  moles  per  liter  in  the  crystal,  and  assuming  an  extinction 
coefficient  of  the  order  of  104  cm-1  M-1  one  can  calculate  a  maximum  path  length  limitation  for 
a  maximum  practical  absorbance  value  of  3  to  be  about  10-4  cm.  Crystals  of  these  oligomer 
duplexes  are  "blocky"  in  nature,  and  direct  absorption  experiments  with  micron  sized  crystals 
would  be  very  difficult  at  best  In  order  to  avoid  the  thin  crystal  requirement  of  direct  absorption 
measurements,  we  have  measured  the  polarized  reflection  spectra  from  the  be  crystal  face  of 
relatively  large  single  crystals  (-200  pi).  The  corresponding  absorption  curves  are  subsequently 
obtained  by  Kramers-Kronig  analysis  of  the  reflection  data.  The  instrument  used  and  the 
transformation  procedure  have  been  described  in  other  publications.4  A  15x  Cassegrainian 
reflective  microscope  objective  is  used  to  focus  the  incident  radiation  onto  the  crystal  face.  The 
incident  beam  is  sharply  tailored  to  minimize  depolarization  effects  and  to  reduce  the  convergent 
nature  of  the  focused  light  to  an  average  angle  of  1°  from  normal  incidence.  A  100  x  200  micron 
image  of  the  monochromator  slit  is  formed  on  the  crystal  surface.  The  dimensions  of  the  crystal 
face  must  be  somewhat  larger  than  that  image. 


For  the  spectral  measurements  the  crystals  were  mounted  on  a  small  pedestal  inside  a  delrin 
capsule  which  incorporated  a  demountable  window.  A  small  drop  of  mother  liquor  (50:50 
waten2-methyl-2,4-pentadiol)  was  used  to  adhere  the  crystal  to  the  pedestal.  Additional  drops 
of  mother  liquor  were  placed  inside  the  capsule  to  maintain  a  saturated  atmosphere  around  the 
crystal.  The  capsule  was  attached  to  a  goniometer  head  and  the  crystal  face  brought  normal  to 
the  incident  light  beam.  The  whole  assembly  was  rotated  about  an  axis  colinear  with  the 
incident  light  so  as  to  align  the  chosen  crystal  axis  with  the  polarization  of  the  light  beam. 
Finally,  the  window  was  removed,  the  alignment  readjusted  as  needed  and  the  reflection  spec¬ 
trum  scanned  from  370  to  190  nm.  The  crystal  was  subjected  to  the  "drying"  effects  of  the 
laboratory  air  for  about  40  minutes.  No  degradation  was  observed,  however,  during  the  course  of 
the  40  minutes  nor  during  a  similar  measurement  period  on  a  subsequent  day  about  two  months 
later. 

The  reflection  spectra  of  the  be  face  are  shown  in  Fig.  4.  For  such  small  crystals  (-0.1  mm 
thick)  mounted  with  the  backside  in  contact  with  a  film  of  mother  liquor,  one  would  expect  there 
would  be  "extra"  reflections  that  add  to  the  signal  throughout  the  transparent  region.  The  extra 
reflections  arise  from  the  crystal/mother  liquor  boundary  and  the  mother  liquor/delrin  (pedestal) 
boundary.  These  additional  reflections  are  of  small  magnitude  relative  to  the  reflection  from  the 
air/crystal  boundary  (front  surface)  because  the  indices  of  refraction  of  crystal,  mother  liquor  and 
delrin  are  similar  in  magnitude.  However,  upon  the  onset  of  absorption  in  the  crystal,  virtually 
no  light  manages  to  penetrate  as  far  as  the  back  surface  of  the  crystal,  and  the  extra  contributions 
vanish.  The  dashed  extensions  of  the  reflectivity  curves  in  Fig.  4  represent  our  attempt  to  remove 
these  unwanted  reflections.  The  onset  of  absorption  in  these  crystals  is  near  310  nm  and  the 
extensions  were  drawn  to  smoothly  connect  with  the  curves  at  the  point  where  only  the  front  cry¬ 
stal  surface  contributes  to  the  signal.  The  use  of  the  correct  reflectivity  in  the  transparent  region 
leading  up  to  the  onset  of  absorption  is  important  in  the  Kramers-Kronig  procedure,  for  it  is  the 
calculated  outcome  in  this  region  that  is  used  as  a  criterion  for  the  adequacy  of  the  transforma¬ 
tion.  The  transformation  involves  numerical  integrations  over  the  entire  frequency  domain 


(0— >°°),  whereas  measurements  are  made  over  a  "somewhat",  narrower  region  (370-190  nm). 
The  usual  procedure  is  first  to  extend  the  measured  spectrum  with  trial  reflectivities  on  the  red 
side  (to  zero  frequency)  and  in  the  vacuum  UV  (to  very  large  but  finite  frequency).4  Then  these 
extensions  are  adjusted  in  a  systematic  way  until  zero  absorpdon  is  calculated  throughout  the 
known  transparent  region.  The  result  is  not  unique  and  depends  on  the  particular  shape  assumed 
for  the  adjacent  region  in  the  vacuum  UV.  By  employing  a  variety  of  vacuum  UV  extensions  a 
family  of  absorption  curves  can  be  obtained  that  mutually  satisfy  the  transparent  region  con¬ 
straint.  The  derived  absorption  spectra  are  shown  in  Fig.  5.  The  range  of  values  for  the  computed 
extinction  coefficients  that  arises  from  the  various  vacuum  UV  trial  functions  are  indicated  by 
the  widening  fans.  The  curves  show  broad,  diffuse  bands  with  no  apparent  structure. 

The  removal  of  the  unwanted  extra  reflection  in  the  transparent  region  (the  dashed  exten¬ 
sions  of  Fig.  4)  is  warranted  for  two  reasons.  First,  different  crystals  show  different  magnitudes 
of  extra  reflection,  however  in  all  cases  the  effect  disappears  at  310  nm.  Second,  the  raw  data, 
i.e.  the  uncorrected  reflection  curve,  cannot  be  satisfactorily  transformed.  The  step  in  the 
reflection  at  320  nm  always  results  in  an  absorption  spectrum  with  a  negative  dip  at  370  nm  and 
a  positive  peak  at  323  nm.  Our  experience  over  some  years  with  the  Kramers-Kronig  transforma¬ 
tion  of  reflection  data  leads  us  to  assert  that  a  satisfactory  transform  cannot  be  obtained  from  a 
reflection  spectrum  that  has  an  incorrect  shape.  We  conclude  that  the  low  energy  extensions 
shown  in  Fig.  4  are  appropriate. 

The  most  likely  source  of  error  in  the  computed  extinction  coefficients,  in  our  estimate, 
arises  from  imperfections  on  the  crystal  surface.  Surface  imperfections  tend  to  scatter  the 
incident  radiation  and  uniformly  lower  the  reflected  signal.  Again,  our  experience  with  the 
Kramers-Kronig  analysis  of  reflection  data,  suggests  that  a  given  uniform  percentage  lowering  of 
the  reflection  produces  an  equivalent  percentage  error  in  the  computed  absorption  spectrum.  The 
correct  shape  is,  however,  maintained.  The  crystal  surfaces  used  in  this  study  showed  regions 
that  were  smooth  and  shiny  under  microscopic  examination.  Although  these  regions  ought  to 
have  yielded  proper  reflectivity,  no  independent  check  of  the  accuracy  of  the  reflection 
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coefficients  was  available. 


Interpretation 

For  our  purposes  here  we  will  assume  that  the  purine  and  pyrimidine  bases  are  the  only 
chromophores  that  directly  contribute  to  the  absorption  at  wavelengths  greater  than  190  nm.  The 
deoxyribose  moieties,  the  phosphate  linkages,  as  well  as  the  water  molecules  present  in  the 
structure  are  not  expected  to  absorb  at  wavelengths  longer  than  about  180  nm.  The  excited  elec¬ 
tronic  states  of  the  individual  bases  are  now  known  fairly  well  (vide  infra)  and  conform  with  the 
expectations  of  basic  iz  electron  theory.  The  strongest  bands  (£ave  =»  204cm-1  M-1)  are  assigned 
to  tc— »tc*  transitions  and  are  all  polarized  in  the  molecular  plane,  n— »tt*  bands  are  expected  to  be 
polarized  normal  to  the  planes  of  the  bases,  but  little  evidence  for  this  type  of  transition  has  been 
found.  Any  such  bands,  if  present,  ought  to  be  weak  and  presumably  are  lost  beneath  the  much 
stronger  in-plane  polarized  bands. 

On  the  other  hand  the  model  for  an  electronic  band  as  being  represented  by  a  single  transi¬ 
tion  moment  direction  is  an  approximation.  Through  vibrational  distortions  the  individual 
vibronic  transitions  may  appear  with  different  polarizations.  Those  that  remain  polarized  in¬ 
plane  will  appear  with  slightly  varying  in-plane  directions;  they  should,  however,  cluster  around 
some  dominant  direction.  Some  vibrational  distortions  lead  to  "forbidden  components"  polar¬ 
ized  normal  to  the  molecular  plane.  The  result  is  that  even  a  nominally  in-plane  polarized  7C — >>7C* 
transition  can  be  expected  to  show  some  out-of-plane  transition  probability.  Transition  moment 
data  acquired  from  spectral  measurements  on  single  crystals  of  individual  bases  suggest  that  3- 
5%  of  a  nominally  in-plane  polarized  intensity  can  be  expected  to  appear  polarized  normal  to  the 
molecular  plane.4""6 

Returning  to  the  crystal  absorption  spectra  for  the  37,000  cm-1  (270  nm)  band  we  see  the  b 
axis  spectrum  (i.e.  close  to  parallel  to  the  planes  of  the  bases)  shows  an  absorptivity  of  9,100 
cm-1  M-1  whereas  the  lie  spectrum  exhibits  a  corresponding  value  of  about  500.  The  component 
oscillator  strengths  evaluated  with  the  formula  f  =  4.32  x  10-5  J  e  d  v  where  e  is  the  decadic 
molar  extinction  coefficient  and  v  is  the  frequency  in  cm-1  are  0.26  and  0.017  for  the  b-axis  and 
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'  c-axis,  respectively.  The  oscillator  strength  appearing  in  the  c.  direction  is  about  three  times 
larger  than  would  be  anticipated  if  it  arose  solely  from  the  8°  cant  of  in-plane  transition  moments 
of  the  bases.  The  difference  along  c  corresponds  to  a  few  percent  of  the  in-plane  intensity  and, 
as  discussed  above,  likely  corresponds  to  the  out-of-plane  components  of  the  nominally  in-plane 
polarized  7t-rc*  transitions  rather  than  perpendicularly  polarized  nrc*  bands.  In  any  event,  the 
dichroic  ratio  is  fairly  large,  (ff/fc)  =  15.  Given  all  the  above  considerations  we  see  that  the 
experimental  spectra  conform  to  expectations  in  their  general  appearance. 

Model  Calculations 

Here  we  will  calculate  the  spectrum  of  the  duplex  from  known  monomer  transition  proper¬ 
ties  and  the  known  relative  positions  and  orientations  of  the  bases  in  the  crystal.  The  simplest 
approach  involves  treating  the  duplex  as  an  "oriented  gas"  of  twelve  bases  each  of  which  will  be 
carrying  its  personal  load  of  electronic  transitions.  Each  transition  moment  is  projected  onto  the 
a,  b  and  c  crystal  axes  and  the  component  oscillator  strength  then  computed.  For  comparison  to 
experimental  crystal  spectra  the  individual  components  bands  are  convened  to  Gaussian  shaped 
bands  using  a  half-width  that  is  common  to  various  crystal  spectra  of  the  individual  bases  (4,500 
cm-1)  and  then  summed  to  yield  the  final  computed  absorption  spectra  For  comparison  to  ran¬ 
dom  solution  spectra  the  transition  moment  data  are  used  directly  to  compute  a  summed  Gaus¬ 
sian  spectrum. 

The  optical  effects  of  base-base  interactions  can  be  handled  most  simply  using  an  Indepen¬ 
dent  Systems  or  Exciton  Model  approach  and  employing  the  point  dipole  approximation  to 
evaluate  interaction  terms.  Conventional  wisdom  suggests  that  the  point  dipole  approximation  is 
not  adequate  for  the  nearest  neighbor  interactions  along  a  polymer  chain  and  that  distributed 
monopolcs  (or  higher  multipoles)  ought  to  be  used.7  The  trouble  in  this  regard  is  that  the  mono¬ 
pole  distribution  for  each  transition  must  be  obtained  from  LCAO  wave  functions  of  dubious 
trustworthiness.  At  the  present  time,  there  is  no  satisfactory  assignment  of  the  experimentally 
observed  monomer  transitions  according  to  any  of  the  various  LCAO  quantum  mechanical 

it 

results.  Treating  the  interactions  in  the  poinr  dipole  approximation  using  all  experimentally 
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determined  transition  moment  directions  seems  to  be  an  appropriate  compromise  to  take  at  the 
present  time.  The  excited  state  eigenfunctions  of  the  duplex  are  obtained  as  linear  combinations 
of  all  monomer  excited  states  by  diagonalizing  a  Hamiltonian  matrix  whose  off  diagonal  ele¬ 
ments  are  mathematically  equivalent  to  the  electrostatic  interaction  of  the  two  transition  moment 
dipoles  from  the  ground  state  to  the  states  involved  in  that  matrix  element.  These  terms  are 
given  by  the  usual  expression 


hia,j{5  —  (1 
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where  a  and  {3  are  the  transitions  on  sites  i  and  j  respectively.  The  d’s  are  transition  moment 
vectors  and  is  the  vector  between  the  respective  base  centers  (i— >j).  The  Kronecker  delta,  5^, 
eliminates  terms  between  transitions  at  the  same  site,  and  diagonal  terms,  hia>ia  are  just  the 
unperturbed  transition  energies.  New  (perturbed)  transition  moment  vectors  are  evaluated  from 
the  resultant  eigenfunctions  and  projected  individually  onto  the  a,  b  and  c  crystal  axes  at  their 
computed  new  energies.  The  results  are  again  Gaussianized  and  finally  presented  as  spectral 
curves  appropriate  for  a,  b ,  c  crystal  axis  spectra.  There  are  74  transitions  covering  a  range  from 
33,000  to  74,000  cm-1  (300-135  nm)  to  be  included  in  the  twelve  bases  of  the  duplex  so  that  the 
matrix  to  be  diagonalized  is  74  x  74.  Gaussianizing  the  output  seems  an  appropriate  tack  to  take 
in  presenting  and  comparing  these  calculated  results  owing  to  the  complexity  of  the  assembly 
and  the  large  number  of  overlapping  transitions  that  must  be  dealt  with. 

Transition  moment  data  for  the  individual  transitions  of  the  four  bases  have  been  obtained 
mostly  from  single  crystal  spectra  of  the  bases  (or  various  methylated  forms).4,5,9-12  We  sum¬ 
marize  these  data  in  Table  II.  In  certain  higher  energy  cases  an  ambiguity  remains  in  the  crystal 
result  and  the  choices  excluded  in  the  calculations  reported  here  are  indicated  by  parentheses. 
For  this  reason,  the  calculated  spectra  in  the  210-190  nm  region  are  presented  provisionally. 
Band  contours  in  this  wavelength  region  do  depend  on  the  various  choices  possible,  particularly 
for  the  fourth  cytosine  transition.  In  order  to  proceed  in  a  consistent  manner  we  have  also  used 
transition  energies  taken  from  crystal  spectra.  The  transition  energy  differences  between  solution 
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and  crystal  environments  is  usually  not  large  in  any  event.  With  the  possible  exception  of  the 
46,000  cm-1  (217  nm)  band  of  the  uracil  chromophore,  no  perpendicularly  polarized  transidon 
has  been  observed  with  any  certainty.  Most  of  the  perpendicularly  polarized  intensity  observed 
in  crystal  spectra  is  more  likely  to  be  simply  "forbidden''  components  of  the  nominally  in-plane 
polarized  transitions. 

The  results  of  the  calculations  for  the  a ,  b  and  c  axes  are  presented  in  Fig.  6  along  with  the 
experimental  spectra  for  comparison.  Both  schemes  (i.e.  oriented  gas  and  with  interactions) 
predict  noticeable  differences  in  the  curves  along  the  a  and  b  axes  particularly  in  the  210-190  nm 
range.  The  first  band  maximum  appears  at  265  nm  along  b  and  at  261  nm  along  a  in  the  oriented 
gas  analysis.  After  base-base  interactions  are  included  there  is  a  slight  blue  shift  of  the  peak 
position  of  the  a  axis  spectrum  (to  259  nm)  and  an  even  smaller  red  shift  along  the  6-axis.  These 
changes  ought  not  be  considered  only  as  shifts  of  energy  levels  but  also  as  resulting  from  the 
redistribution  of  intensity  among  the  multitude  of  closely  spaced  (overlapped)  transitions. 

The  theory  of  the  hypochromic  effect  in  stacked  polynucleotides  was  worked  out  by 
Tinoco13,14  and  by  Rhodes.13  The  in-plane  transition  moments  of  an  approximately  card-pack 
stack  of  bases  mix  under  their  mutual  interactions  so  as  to  shift  intensity  from  lower  energy 
states  to  higher  energy  states.  Transition  moments  polarized  normal  to  the  planar  bases  show  the 
opposite  effect  (hyperchromism)  owing  to  their  approximate  head-to-tail  arrangement.  Model 
calculations  in  the  point  dipole  approximation  presented  by  both  Tinoco  and  Rhodes  in  the  early 
work  were  done  with  assumed,  artificial  transition  moment  data,  and  the  results  were  simply 
shown  to  be  consistent  with  the  observed  sign  and  general  magnitude  of  the  observed  effects. 
More  recent  work  by  Rizzo  and  Schellman,7  for  example,  use  a  mix  of  theoretical  and  empirical 
parameters  for  the  monomer  transition  moment  data  set.  To  our  knowledge,  however,  there  has 
been  no  hypochromism  calculation  at  any  level  of  approximation  in  which  only  experimental 
transition  moment  data  is  used  and  that  brings  out  the  spatial  anisotropy  of  the  duplex.  The  cal¬ 
culations  reported  here  pertain  to  incident  light  polarized  along  the  a,  b  and  c  crystal  axes  of  the 
duplex,  and  since  the  spectral  properties  are  not  the  same  in  these  different  directions,  the  results 
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contain  more  verifiable  information  than  the  outcome  of  isotrop.ic  calculations.  Of  course,  for 
comparison  to  isotropic  solution  spectra  the  calculated  results  can  be  randomized  by  summing 
the  a,  6  and  c  axis  spectra  and  dividing  by  3. 

The  calculated  curves  shown  in  Fig.  6  give  18%,  16%  and  6%  reductions  in  oscillator 
strength  along  the  n,  b  and  c  crystal  axes  respectively  for  the  260  nm  region.  Only  the  experi¬ 
mental  6-axis  result  is  available  for  comparison  to  its  corresponding  calculated  value,  and  here 
the  observed  loss  of  absorption  of  45%  corresponds  to  about  2.8x  more  reduction  than  that  cal¬ 
culated.  Along  both  a  and  b  axes  there  is  calculated  hypochromicity  throughout  the  low  energy 
region,  while  the  higher  energy  region  (220-180  nm)  shows  mixed  hypo/hyperchromicity. 
Although  not  apparent  in  Fig.  6,  the  very  weak  absorption  intensity  polarized  parallel  to  the  helix 
axis  shows  a  slight  decrease  in  f  number  for  the  low  energy  region  arising  from  base-base 
interactions.  This  result  occurs  since  the  parallel  spectrum  arises  simply  as  the  components  of 
the  in-plane  polarized  transitions  that  are  undergoing  hypochromic  changes. 

In  order  to  investigate  the  difference  between  the  theoretical  and  experimental  results  for 
the  6-axis  hypochromism  we  compare  in  Fig.  7  the  randomized  crystal  spectrum  to  the  aqueous 
solution  spectrum  of  an  equivalent  mix  of  the  constituent  deoxynucleosides.  Although  the  curve 
derived  for  the  randomized  crystal  spectrum  using  the  oriented  gas  approach  differs  in  shape 
somewhat  from  that  of  the  nucleoside  mix,  the  oscillator  strengths  differ  by  only  7%  (f^o  =  0-32 
for  the  mix  and  0.30  for  randomized  oriented  gas).  The  comparison  is  judged  reasonable  consid¬ 
ering  the  transition  moment  data  are  taken  from  various  related  purine  and  pyrimidine  bases  and 
not  the  specific  deoxynucleosides  present.  Since  there  are  twelve  base  pairs  per  helical  turn  one 
might  guess  the  a  and  6  axes  to  show  similar  spectra,  and  this  supposition  is  supported  by  the 
calculations  which  yield  only  -9%  more  intensity  along  a  than  along  6  and  similar  values  for  the 
calculated  hypochromism.  Assuming  such  rough  isotropy  is  the  case,  we  can  (in  the  absence  of 
experimental  a  axis  data)  construct  an  approximate  randomized  crystal  spectrum  by  taking  twice 
the  6  axis  spectrum,  adding  the  c  axis  spectrum  and  dividing  by  3.  The  result  is  a  curve  where 
=  270  nm,  =  6,200  cm*1  M"1  and  f  (isotropic)  =  0.18.  The  last  two  numbers  may  be 
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-5%  low  judging  from  the  9%  more  absorption  predicted  along  the  a  axis  compared  to  the  b- 
axis.  The  isotropic  hypochromism  [H  =  (1  -  of  the  long  wavelength  band  upon 

including  base-base  interactions  is  0.17,  whereas  the  approximate  isotropic  spectrum  derived 
from  the  crystal  data  gives  a  value  of  H  =  0.45  (in  comparison  to  the  mix  of  deoxynucleosides). 

As  mentioned  earlier,  surface  imperfections  may  contribute  to  the  measured  hypochrom¬ 
ism.  However,  another  source  of  the  above  difference  in  H  may  arise  in  chain  length  considera¬ 
tions.  It  must  be  recognized  that  the  stacking  arrangement  of  duplexes  in  the  crystal  generates 
essentially  continuous  double  helical  chains,  whereas  the  calculations  include  only  the  interac¬ 
tions  within  a  single  duplex.  A  more  appropriate  comparison  for  these  calculated  results  would 
seem,  therefore,  to  be  with  solution  spectra  of  the  duplex  itself.  Fig.  8  shows  the  aqueous  solu¬ 
tion  spectra  of  the  melted  form  (80°C)  and  both  B  and  Z  double  helical  forms.  Relative  to  the 
melted  curve  the  260  nm  band  exhibits  a  value  of  H  =  0.19  for  the  Z-form,  and  this  value  com¬ 
pares  reasonably  with  the  calculated  result  of  0.15.  Such  agreement  tends  to  support  the  asser¬ 
tions  that  (i)  the  crystal  spectra  pertain  to  essentially  infinite  helices  and  (ii)  the  calculation  of  the 
absorptive  properties  of  polynucleotides  may  be  adequately  approached  within  the  framework  of 
the  point  dipole  approximation.  The  calculations  we  have  done  lead  to  no  hyperchromicity  in  the 
long  wavelength  region  (-295  nm),  whereas  such  is  observed.  Perpendicularly  polarized,  mt* 
intensity  has  been  presumed  to  provide  the  mechanism  for  this  effect.  However,  our  calculations 
included  no  such  transition  in  the  long  wavelength  region.  The  weakness  of  the  experimental 
absorption  polarized  parallel  to  the  helix  axis  (c  crystal  axis)  rules  out  the  mt*  mechanism  for 
the  hyperchromicity  in  the  long  wavelength  region. 

SUMMARY 

Absorption  spectra  of  the  be  crystal  face  of  single  crystals  of  the  Z-form  duplex 
d(m3CGUAm3CG)-d(m3CGUAm3CG)  have  been  obtained  from  polarized  reflecdon  data.  The 
crystal  structure  is  such  that  the  results  arc  appropriate  for  an  infinite,  Z-form  double  helix  with  a 
repeaung  unit  of  two  six  base-pair  duplexes.  The  c  axis  spectrum  is  parallel  to  the  helix  axis  and 
shows  very  small  absorption  intensity,  whereas  the  b-axis  spectrum  is  perpendicular  to  the  helix 
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axis  and  shows  the  7t7t*  base  transitions  strongly. 

Calculations  based  on  both  an  oriented  gas  model  and  a  point  dipole  interaction  model  sug¬ 
gest  that  the  ab  crystal  face  (i.e.  perpendicular  to  the  helix  axis)  will  show  significant  differences 
between  the  a  and  b  axis  spectra.  The  results  for  the  dipole  interaction  model  lead  to  a  hypo- 
chromism  of  0.15  compared  to  an  experimental  value  of  0.45.  The  disparity  is  most  likely  attri¬ 
butable  to  the  chain  length  differences  between  the  crystal  (infinite)  and  duplex  (six-mer). 
Efforts  at  preparing  crystals  with  (ab)  faces  large  enough  for  measurement  and  extending  the 
spectra  deeper  into  the  vacuum  LTV  are  presently  underway. 

*This  work  was  supported  by  grants  from  the  NIH  (#GM38575  to  L.B.C.),  NIEHS  (#E500210  to 
P.S.H.)  and  ONR  (#N00014-88-K-0388  to  P.S.H.). 
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TABLE  I 


Coordinates  of  Molecule  Centers  and  Direction  Cosines  of 
the  Reference  Axes  and  Normals  to  the  Planes 


Center1 

Reference  Axisb 

Normal1* 

a 

b 

c 

a 

b 

c 

a 

b 

c 

Cl 

0.6847 

0.3769 

-0.0765 

0.5728 

-0.8139 

-0.0973 

0.1486 

-0.0136 

0.9888 

G2 

0.4368 

0.4481 

0.0179 

0.0619 

-0.9948 

0.0805 

0.0540 

0.0839 

0.9950 

U3 

0.5540 

0.4986 

0.0938 

-0.4433 

-0.8963 

0.0147 

0.0389 

-0.0029 

•  0.9992 

A4 

0.5245 

0.6505 

0.1771 

-0.8254 

-0.5594 

-0.0753 

-0.0862 

-0.0069 

0.9963 

C5 

0.6102 

0.6101 

0.2581 

-0.9904 

-0.1261 

-0.0572 

-0.0482 

-0.0734 

0.9961 

G6 

0.8292 

0.6962 

0.3429 

-0.8760 

0.4810 

-0.0349 

-0.0538 

-0.0256 

0.9982 

C7 

0.5524 

0.5575 

0.3310 

-0.1438 

0.9868 

0.0741 

0.0899 

0.0876 

-0.9921 

G8 

0.4574 

0.4126 

0.2438 

-0.6699 

0.7423 

-0.0168 

0.0445 

0.0176 

-0.9989 

U9 

0.6065 

0.4390 

0.1669 

-0.9480 

0.3153 

0.0427 

j 

-0.0110 

0.1017 

-0.9948 

A10 

0.7570 

0.3220 

0.0822 

-0.9585 

-0.2635 

0.1086 

-0.0912 

-0.0775 

-0.9928 

Cll 

0.7826 

0.3779 

-0.0024 

-0.7854 

-0.6083 

0.1146 

-0.1810 

0.0486 

-0.9823 

G12 

1.0520 

0.3958 

-0.0886 

-0.3025 

-0.9530 

0.0133 

0.0111 

-0.0175 

-0.9998 

^e  centers  of  the  cytosine  and  uracil  residues  were  taken  as  the  midpoint  between  the  Nt  and 
C4  line  segment  while  the  midpoint  of  the  C4-C5  bond  was  taken  as  the  center  of  the  guanine 
and  adenine  residues.  These  coordinates  are  given  as  fractions  of  the  unit  ceil  dimensions. 


b  The  reference  axis  is  the  C4-C5  line  in  the  purines  and  the  N1-C4  line  in  the  pyrimidines.  The 
normal  was  constructed  by  taking  the  vector  product  of  the  reference  axis  vector  with  the  vector 
-  directed  from  C4  to  Cj  in  the  purines  and  Ni  to  C2  in  the  pyrimidines. 


-77- 


TABLE  II 

Monomer  Transition  Moment  Data 


va  (cm-1) 

X  (nm) 

f 

(A) 

0b 

Ref 

Adenine 

1 

38,500 

260 

0.11 

0.513 

83 

16 

2 

37,500 

267 

0.20 

0.701 

25 

3 

47,000 

213 

0.25 

0.700 

-45 

4 

49,000 

204 

0.11 

0.455 

15 

5 

53,700 

186 

0.30 

0.718 

72 

12 

6 

56,000 

179 

0.10 

0.406  •— 

~45 - 

7 

62,500 

160 

0.23 

0.582 

6 

8 

69,900 

143 

0.10 

0.363 

-45 

Cytosine 

1 

36,900 

271 

0.14 

0.591 

6 

2 

43,000 

233 

0.03 

0.254 

-35 

3 

44,600 

224 

0.13 

0.518 

76 

4 

4 

50,500 

198 

0.36 

0.811 

86  (-26)c 

5 

62,000 

161 

0.15 

0.472 

0(60) 

6 

66,000 

152 

0.20 

.  0.529 

60(0) 

Guanine 

1 

36,000 

272 

0.16 

0.640 

-4 

2 

39,400 

254 

0.25 

0.765 

-75 

3 

49,000 

206 

0.41 

0.878 

-71 

10 

4 

53,000 

187 

0.48 

0.914 

41 

5 

63,000 

159 

0.20 

0.541 

0(40) 

6 

69,000 

145 

0.20 

0.517 

40  (0) 

Uracil 

1 

38,500 

260 

0.19 

0.674 

-9 

9 

2 

47,000 

213 

0.26 

0.714 

-53(+59)d 

5 

3 

46,000 

217 

0.007 

0.118 

1 

4 

55,000 

182 

0.31 

0.721 

-33(+39°) 

17 

5 

60,000 

167 

0.17 

0.511 

-40(+47) 

•Taken  as  peak  positions  in  the  crystal  spectra  of  individual  bases. 


bIn-plane  transition  moment  directions  using  the  Tinoco-DeVoe  convention.  For  pyrimidines  9 
is  positive  toward  N3  from  the  N1-C4  reference  axis,  and  for  the  purines  0  is  positive  toward  N3 
from  the  C4-C5  reference  axis.  Values  in  parentheses  are  alternate  assignments  owing  to  ambi¬ 
guities  in  the  crystal  spectra  of  the  particular  base. 

®The  ambiguity  in  the  polarization  direction  of  transition  4  of  cytosine  has  not  been  resolved. 
The  +86°  choice  provides  for  a  good  fit  of  both  CD  and  LD  spectra  of  various  G  and  C  contain¬ 
ing  polymers,  while  the  -26°  choice  does  not  lead  to  a  satisfactory  fit.18  On  the  other  hand, 
theoretical  calculations  uniformly  predict  values  in  the  -35°  to  -45°  range. 

dThe  ambiguity  for  band  2  had  not  been  satisfactorily  resolved.  The  -53°  choice  has  been  used 
to  fit  CD  and  LD  spectra  of  polymers  containing  A  and  T.19  On  the  other  hand,  theoretical  calcu¬ 
lations  invariably  yield  values  -50°. 8  Additionally  the  +59°  direction  harmonizes  better  with 
LD  studies  using  stretched  polyvinylalcohol  sheets. 


Figure  Captions 


Fig.  1.  Projection  of  six  duplexes  onto  the  be  face.  The  helix  axes  are  parallel  to  the  c  crystal 
axis.  The  four  unit  cell  sites  l(x,y,z),  2(x+y,  y,  z+y),  3(x,y  +  y,z+y)  and  4(x+y,y+y,z)  are  indi¬ 
cated  in  the  figure.  The  base  numbering  scheme  runs  1  to  6  down  one  chain  of  m5CGUAm5CG 
and  7-12  up  the  complementary  sequence  of  the  duplex.  The  planes  of  the  bases  are  tilted  by  an 
average  8°  from  the  ab  crystal  plane. 

Fig.  2.  Projection  onto  the  be  crystal  plane  using  the  van  der  Waal’s  surfaces  of  the  atoms.  The 
shaded  duplexes  are  at  the  backside  of  the  unit  cell,  whereas  the  unshaded  duplex  runs  through 
the  center  of  the  unit  cell.  The  continuous  nature  of  the  end-to-end  stacking  of  duplexes  is 
brought  out  in  this  figure.  The  solid  lines  connect  the  phosphorus  atoms  of  the  phosphate  back¬ 
bone.  Atoms  designations  are  as  follows:  closed  concentric  circles  =  carbon,  dashed  concentric 
circles  =  nitrogen,  dotted  concentric  circles  =  oxygen,  saw-toothed  concentric  circles  -  phos¬ 
phorus. 

Fig.  3.  Crystal  morphology  of  the  Z-DNA  duplex  of  d(m5CGUAm5CG).  The  be  face  [with 
Miller  indices  (100)]  was  the  face  examined  in  this  study.  The  dimensions  are  about 
0.2  x  0.2  x  0. 1  mm. 

Fig.  4.  Reflection  spectra  for  radiation  polarized  parallel  to  the  c  and  b  crystal  axes.  These  axes 
are  parallel  and  perpendicular  to  the  helix  axis  respectively.  The  dashed  portions  are  the  exten¬ 
sions  used  in  the  transparent  region  for  the  Kramers-Kronig  transformation  procedure.  These 
dashed  extensions  are  drawn  so  as  to  eliminate  "extra"  reflection  contributions  from  the  back 
side  of  the  crystal  in  the  transparent  region  as  discussed  in  the  text. 

Fig.  5.  Absorption  curves  for  the  b  and  c  axis  derived  from  the  reflection  spectra  shown  in  Fig.  4 
by  Kramers-Kronig  analysis.  The  diverging  width  of  the  curves  toward  higher  cm-1  correspond 
to  the  variations  found  for  different  trial  vacuum  UV  reflectivities  used  in  the  transformation 
procedure.  The  "wiggles"  arise  from  instrumental  noise. 
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Fig.  6.  Comparison  of  calculated  and  experimental  crystal  spectra.  (Top)  Calculated  a-axis 
spectrum  with  (dots)  and  without  (solid)  interbase  interactions.  (Middle)  Calculated  6-axis  spec¬ 
trum  with  (dots)  and  without  (solid)  interbase  interactions.  The  experimental  6-axis  curve  from 
Fig.  5  is  shown  for  comparison.  (Bottom)  The  various  curves  are  as  above,  but  here  represent  the 
c-axis.  Note  the  expanded  vertical  scale.  Even  on  this  scale  the  difference  between  the  oriented 
gas  and  interactions  models  is  not  discernible. 

Fig.  7.  Comparison  of  isotropic  or  randomized  spectra.  The  ordinate  scale  is  on  a  per  mole  base 
chromophore  basis.  The  solid  curve  (a)  corresponds  to  the  absorption  spectrum  of  an  equivalent 
mixture  of  the  deoxynucleosides  in  aqueous  solution.  The  dashed  curve  (c)  represents  the  ran¬ 
domized  spectrum  of  the  transition  moment  data  used  for  the  four  types  of  chromophores  in  the 
theoretical  calculations  and  referred  to  as  the  oriented  gas  model  approach  in  the  text.  The  dotted 
curve  is  the  randomized  result  after  base-base  interactions  have  been  included.  The  solid  curve 
was  obtained  from  an  equivalent  mix  of  deoxy nucleotides  in  neutral  aqueous  solution. 

Fig.  8.  Solution  spectra  of  d(m5CGUAm5CG)  in  various  conformations.  The  dotted  curve  is  for 
the  random  coil  or  melted  (80°C)  form.  “The  dashed  curve  represents  a  duplex  in  the  B-  - 
conformation  while  the  dash-dot  curve  corresponds  to  the  Z-form  duplex.  Since  actual  duplex 
concentrations  were  unknown,  the  spectrum  for  the  melted  form  were  scaled  to  match  the  max¬ 
imum  of  the  deoxynucleoside  mix.  The  B-  and  Z-form  spectra  were  scaled  by  the  same  factor. 
Although  there  is  residual  hypochromism  upon  melting  of  long  duplexes,  there  should  be  little 
such  hypochromism  for  hexamers.  The  B-DNA  spectrum  was  recorded  in  0.1  mM  MgCl2  at  pH 
7  and  10°C.  The  B  conformation  was  confirmed  by  CD  measurements.  The  spectrum  of  the 
melted  form  was  obtained  simply  by  raising  the  temperature  to  80°C.  The  Z-conformer  could  not 
be  obtained  in  pure  form  but  only  as  a  mixture  with  the  B-form.  The  difference  spectrum 
between  the  B-form  and  the  mixture  was  used  along  with  known  differences  between  the  B  and 
Z  forms  of  poly(m5CG)  to  construct  the  spectrum  representing  the  pure  Z-form.  The  values  for 
the  hypochromism  of  the  260  run  band  obtained  from  these  curves  are  0.16  and  0.19  for  the  B- 
and  Z-forms  respectively.  "  - 
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